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We address three problems faced by effective interfacial Hamiltonian models of wetting based on a 
single collective coordinate £(y) representing the position of the unbinding fluid interface. Problems 
(PI) and (P2) refer to the predictions of non-universality at the upper critical dimension d = 3 at 
critical and complete wetting respectively which are not borne out by Ising model simulation studies. 
(P3) relates to mean-field correlation function structure in the underlying continuum Landau model. 
Building on earlier work by Parry and Boulter we investigate the hypothesis that these concerns arise 
due to the coupling of order parameter fluctuations near the unbinding interface and wall. For quite 
general choices of collective coordinates Xi(y) we show that arbitrary two-field models H[Xi,X2] 
can recover the required anomalous structure of mean-field correlation functions (P3). To go beyond 
mean- field theory we introduce a set TL of Hamiltonians based on proper collective coordinates s(y) 
near the wall which have both interfacial and spin-like components. We argue that an optimum 
model H[s,£] € Tt, in which the degree of coupling is controlled by an angle like variable S* , best 
describes the non-universality of the Ising model and investigate its critical behaviour. For critical 
wetting the appropriate Ginzburg criterion shows that the true asymptotic critical regime for the 
local susceptibility xi is dramatically reduced consistent with observations of mean-field behaviour 
in simulations (PI). For complete wetting the model yields a precise expression for the temperature 
dependence of the renormalized critical amplitude 8 in good agreement with simulations (P2). We 
highlight the importance of a new wetting parameter which describes the physics that emerges due 
to the coupling effects. 

PACS nos: 68.10m, 64.60Fs, 82.65i 



I. INTRODUCTION 

A. Preliminary remarks 

During the last few decades the theory of continuous 
wetting transitions in inhomogeneous fluids and (simple) 
magnets has been extensively developed JO-pl. A 
striking feature of theory is the prediction that the 
critical exponents characterising the phase transition are 
sensitive to the details of the microscopic interactions 
(as well as fluctuation effects). Consider, for example, 
systems with long-range fluid- fluid and wall- fluid forces. 
Above the upper critical dimension, where mean-field 
(MF) theory is valid, the values of the critical exponents 
depend on the precise power law of molecular interac- 
tions Below the upper critical dimension, where 
fluctuations are important, one encounters weak and 
strong fluctuation regimes, with true universality only 
characteristic of the latter . For systems with short- 
range forces the upper critical dimension for both critical 
and complete wetting transitions is d = 3 Below 
this dimension we expect the critical exponents to take 
universal non-classical values at the respective transitions 
j|. However at and above the upper critical dimension 
non-universality may appear in one of two ways. Firstly 
for d > 3, corresponding to the MF regime, critical 
exponents may be non-universal if the decay lengths 
of the (exponential) wall-fluid and fluid-fluid forces are 



different jfj). More interesting is the behaviour at the 
upper critical dimension for which calculations based on 
the standard interfacial (or capillary-wave) model pre- 
dict rather dramatic fluctuation-induced non-universality 
However, as is well known these predictions are not 
supported by extensive Monte Carlo (MC) simulations 
of critical wetting in the three dimensional Ising model 
which are instead consistent with MF theory P,^0[. 
Therefore the status of the theory of wetting in systems 
with short-range forces at the upper critical dimension 
is somewhat controversial — and has been so for over a 
decade. 

In recent years concerted effort has been made to 
refine effective Hamiltonian methods in an attempt to 
overcome this (and other) problem(s) of capillary-wave 
theory for three dimensional systems |TT[ ] . In addition 
developments in theory have been complemented by new 
Ising model simulation studies (f2|,[l3| which raise more 
issues and provide an independent test of the new ideas. 
In this paper we provide comprehensive details of a 
theory of criticality and correlation function structure 
at wetting transitions in systems with short-range forces 
based on a two-field effective Hamiltonian H[sJ] which 
accounts for the coupling of order parameter fluctuations 
at the wall and depinning interface 14 1. The model 
may be regarded as an extension of an earlier coupled 
theory due to Parry and Boulter (PB) Jl5|— |l9[] which 
was largely restricted to the complete wetting phase 
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transition. The new theory, described here also allows 
quantitative discussion of coupling effects at critical 
wetting as well as an improved treatment of the crossover 
from complete to critical wetting. Our analysis is 
consistent with the earlier PB theory which emerges 
naturally as a high temperature limit appropriate to 
modelling complete wetting. We shall argue that the 
coupled theory successfully explains three problems of 
(three dimensional) wetting theory including the critical 
wetting controversy mentioned above. 

As noted above our analysis is a further development 
of the PB theory which is itself a generalisation of a 
model and method introduced by Fisher and Jin (FJ) 
p(j| p5[ . While the FJ model does not allow for coupling 
effects, it doe s g eneralise the standard capillary-wave 
(CW) theory p9] by allowing for a position dependent 
stiffness coefficient which has important consequences. 
We aim to keep our article as self-contained as possible 
and to continue our introduction we review these separate 
developments and focus on a number of important 
issues which provide more detail of the background and 
motivation for the present study. 



B. Microscopic models and the surface phase 
diagram 

A suitable starting point for modelling adsorption 
in continuum fluids (or simple magnets) with short- 
range (contact) forces at a planar wall (situated in the 
z = plane) is the Landau-Ginzburg- Wilson (LGW) 
Hamiltonian fl27u28| 



H L Gw[m] 




dz 



4>(m) 



(1.1) 



based on a local order parameter m(r = (y, z)) which 
we will refer to as the magnetization. We will assume 
throughout this article that the continuum LGW model 
describes the same physics as the Ising model as regards 
wetting properties occurring above the roughening tem- 
perature Tr of the latter (with Tr rs 0.54Tc for a simple 
cubic lattice). The bulk potential function </>(m) yields 
coexistence between bulk phases a and (3 at sub-critical 
temperatures T < Tc and zero bulk field h = 0. We shall 
refer to the bulk magnetizations for the up and down spin 
phases as m a > and mp < respectively including 
cases in which the phase is metastable. The parameters 
h\ and c denote the surface field and enhancement 
respectively while 



mi(y) = m(y,0) 



is the magnetization at the wal l. Often we shall 
abbreviate the surface term in (1.1) as <j)i( m i) = § TO i — 



h\m\. In most of our discussion of correlation functions 
we will let ksT = 1 for convenience but will explicitly 
include it in definitions of wetting parameters and for 
predictions of fluctuation effects. 

The model is expected to show wetting behaviour so 
that for sufficiently large surface field hi > h^(> 0) or 
temperature T > Tw (with Tyy the wetting temperature) 
the wall-/? interface is completely wet by the a phase 
when h = 0~ |^,^]. A section of a possible surface 
phase diagram is sketched in Fig. [l] and shows two types 
of continuous wetting transition. 
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FIG. 1. Section of the surface phase diagram. Along the 
cut in the h — axis the wall-/? interface is completely wet 
by the a phase. In Monte-Carlo simulation studies critical 
wetting has been investigated (for technical reasons) along 
path (ii) , that is with T = Tw and h — > CP . 

At each transition the mean thickness (£) of the ad- 
sorbed a phase diverges leading to large scale fluctuations 
in the position of the a(3 interface characterised by 
perpendicular and parallel correlation lengths £j_ and 
We distinguish between critical wetting transitions which 
occur for T — » T w with fixed surface field (or for fixed 
temperature and hi — > h^~) at h = CP (for example 
routes (i) and (ii) in Fig. ^) and complete wetting in 
which {£} diverges for T > T w (or hi > hf) as h -> 0" 
(shown as route (iii)). If we consider the surface tensions 
involved for the wall-a and wall-/3 interfaces we can define 
a singular contribution /sing for each transition by 



&W/3 = &wa + °ali + /si 



(1.3) 



where a a p is the free fluid interfacial tension between 
the bulk co-existing phases. The wall-/3 interface is 
completely wet by the a phase if / s i n g = corresponding 
to zero contact angle ||. 

For critical and complete wetting transitions occurring 
along routes (i) and (iii) respectively, the main critical 
exponents are defined by 



{I) 



/si 



and 



(1.2) (I) ~ \h\~F a ■ C|| ~ \K 
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(1.4) 



(1.5) 



where r denotes the deviation from the critical wetting 
phase boundary 
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T w -T 



for fixed T 
for fixed h±. 



(1.6) 



A schematic illustration of the various lengthscales in the 
problem is shown in Fig. @. 



wall 

FIG. 2. Schematic picture of the diverging lengthscales at 
a wetting transition. 

In principle the transverse and perpendicular correla- 
tion lengths can be calculated from the pair correlation 
function 

G(r 1; r 2 ) = (m( ri )m(r 2 )) - (m(ri))(m(r 2 )) (1.7) 
or its transverse Fourier transform (with yi 2 = y 2 — yi) 



G(zi,z 2 ;q)= /dyi2e iqyi2 G(n,r 2 ) 



(1.8) 



for particle positions near the a/? interface. 

In effective Hamiltonian theories based on a collective 
coordinate ^(y) representing the local height of the 
interface (at vector position y along the wall) £|| may be 
calculated from the analogous height-height correlation 
function while the root mean square fluctuation is 



a = 



In fact the two lengths are 

not independent of each other and in d = 3 effective 
Hamiltonian theories predict (see [§],§|,f| and also later) 



(1.9) 



(«£l)' = wln 1 + (AM 



with k the inverse (true) correlation length of the bulk a 
phase and A the high momentum cut-off. Here lu is the 
wetting parameter 



k B Tn 2 
4ttS 



a/3 



(1.10) 



with S a( g the stiffness (related to the tension) of the 
free a(3 interface. Importantly d = 3 is the upper 
critical dimension for interfacial roughness (for thermal 
fluctuations) and for d > 3 the interface width ~ £j_ is 
only of order a bulk correlation length. 



A direct analysis of wetting in the LGW model is 
only possible in MF approximation |^,^9| which ignores 
fluctuation effects and is equivalent to an effective value 
of the wetting parameter equal to zero. Minimising (1.1) 
determines the MF free energy 



Fmf = min H LGW [m] 



(1.11) 



and yields a standard Euler-Lagrange equation for the 
MF profile fh(z). For our purposes we note only that for 
c > k |2j| sections of the MF surface phase diagram are 
similar to Fig. [l] with the wetting field and temperature 
determined by the relation 



= cm Q (T, 



MF\ 
W I 



0" 



(1.12) 



The MF values of the critical singularities at critical 
and complete wetting are 



K(t) 



in 



/sing 



and 



K(l)~ln\h\-l ; £||~H- 1/2 ; /, 



sing 



(1.13) 



hln\h\ (1.14) 



respectively. Substitution of the MF critical exponents 
into the hyperscaling relations 2 — a s = (d— l)f|| and 2 — 
a c ° = (d— determines the upper critical dimensions 
as d* = d* D = 3 as quoted earlier It is therefore 

possible that fluctuations alter the value of critical 
exponents and amplitudes in d < 3 but to understand 
this we need to resort to effective Hamiltonian methods. 



C. The capillary- wave model 

The fluctuation theory of wetting is based almost 
entirely on the phenomenological CW model J|] 



H[l] = J dy\h: aP {Vif 



W(£) 



(1.15) 



This has been generally accepted as a suitable coarse- 
grained description of an interface whose fluctuations 
occur on lengthscales much larger than the bulk corre- 
lation length (hence the cut-off restriction A <C K in 
( |l.9[ )). Instead of a microscopic order parameter the 
Hamiltonian is a functional of the collective coordinate 
^(y) mentioned earlier. Long wavelength fluctuations in 
the local position of the interface increase its surface area 
and are resisted by the stiffness term in (1.15). The 
direct or bare interaction of the interface with the wall 
is mediated by the binding potential W(£) constructed 
with reference to MF approaches (see later). For systems 
with short-range forces such as the LGW model the 
form usually assumed for modelling continuous wetting 
transitions is ]7|p6|] 



W{£) = h£ + 2 K m a Te- Ke + be' 



(1.16) 
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together with a hard wall contribution which restricts 
configurations to £(y) > 0. Here h — (rnp — m a )h is 
a measure of the deviation from two-phase coexistence 
while & is a positive constant near Tyy (which is required 
to ensure the stability of the transition). The deviation 
from the critical wetting phase boundary, r, is given 
explicitly as 



hi 



(1.17) 



and implies that the overall sign of the leading exponen- 
tial is positive for T > T^ F . Note that by minimising the 
binding potential we recover the MF expression s for the 
mean film thickness and free energy (1.13, 1.14) quoted 
earlier. Also note that for complete wetting only the first 
exponential term is required to model the transition. 

If we assume that the fluctuations 6£(y) — £(y) — I 
(where (£) = £ in MF theory) are small (which is valid for 
d > 3) then it is straightforward to calculate the structure 
factor 



S(q) = J d yi2e 



«q-yi2 



(S£( yi )6£(y 2 )) 



(1.18) 



which has a simple Lorentzian form 

k B T 



S(q) 



W" + q^ a p 



(1.19) 



reminiscent of Ornstcin-Zcrnikc theory. Here W" 



. which allows us to identify the transverse cor- 
relation £jj = 30 1 . If we assume that small 
fluctuations S£ correspond to magnetization fluctuations 
which translate the MF profile m(z) we are led to the 
prediction |31j 



G(zi,z 2 ; q) « m'(zi)m'(z 2 )S{q) 



(1.20) 



for the MF correlation function for positions z\ , z 2 close 
to the af3 interface. Importantly this agrees with the 
explicit solution of the Ornstein-Zernike equation for the 
LGW model (in MF approximation) which reads |3^ 1 
(setting k B T = l) 



f_cP_ 
\ dz\ 



4>"{fn(zi)) + q 2 ^G(zi 1 z 2 - 1 q) = 8(zi - z 2 ) 

(1.21) 



As we shall see in the next sections MF identifications 



such as (1.20) can be made precise for arbitrary zi,z 2 
using generalized effective Hamiltonian and collective 
coordinate theories. 

For d < 3, renormalization group (RG) and transfer 
matri x studies of the CW model with binding potential 
( 1.16[ ) predict universal critical behaviour for both tran- 
sitions ]l|P,p3[ . For these dimensions the precise form 
of W{£) is not essential because the critical properties 



are determined by fixed points of the RG transformation. 
Thus in d = 2 approximating W by a square well recovers 
the same critical exponents for critical wetting as those 
found in the exact Ising model calculation (34j. However 
for d > 3 the precise structure of W{£) is essential 
because there is no non-trivial fixed point Hamiltonian 

At the upper critical dimension d — 3, linear RG 
calculations predict that the critical singularities are non- 
universal depending on the wetting parameter us [j^). 
For critical wetting the results are particularly dramatic, 
for example the correlation length exponent is given by 




for < uj < | 
for | < uj < 2 
for uj > 2 



(1.22) 



where the last regime corresponds to an exponentially 
fast divergence. Note also that for u> < 2 the wetting 
temperature is unaltered so that r = still denotes 
the critical wetting phase boundary. For complete 
wetting the predictions are less dramatic and critical 
exponents retain their MF values |l. Nevertheless 
critical amplitudes are renormalized and of particular 
interest is the dimensionless adsorption amplitude 



lim 



\n\h\- 



which is given by 



1 + - 

1 T 2 



2w 



for < uj < 2 
for uj > 2 



(1.23) 



(1.24) 



While these predictions are made using an approxi- 
mate linear RG scheme which cannot handle the hard- 
wall contribution to W(£) precisely, there are good 
reasons for believing that the results are exact at least 
in the important regime uj < 2. In particular they are 
supported by MC simulations of a discretized CW model 
|35| and also by numerical studies of a non-linear RG 
analysis which does allow for the hard- wall term 

The value of the wetting parameter is not determined 
by the CW theory and must be regarded as input into 
the model. During the time since these predictions 
were first made the temperature dependence ui(T) for 
the simple cubic Ising model has been studied in some 
detail |36|. Near Tr, uj « 0.5 but its value rises sharply 
on increasing T and is very close to w « 0.8 for all 
temperatures in the range 0.6Tc~Tc- On approaching 
the bulk critical point, hyperuniversality implies that the 
wetting parameter tends to a universal value 



lim 

t->t; 



i LOG 



(1.25) 



estimated as ujc ~ 0.775. Thus for the Ising model 
RG calculations based on the CW model predict the 
numerical values 
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3.7 
1.4 



CW theory 



(1.26) 



(assuming 0.6Tq < j 1 < j 1 ^ significantly different from 
the MF results (G| pi). 



D. The Fisher- Jin model and position dependent 
stiffness 



In a series of important papers 20 - 2% FJ reassessed 
the status of the CW model and suggested specific 
modifications that have some interesting consequences 
in = 3. They were principally motivated by one 
of the central problems of wetting theory (see (PI) 
below) and forwarded a novel explanation. In retrospect 
however there are further problems that need to be 
addressed that the FJ model cannot provide answers 
for and it is our view that the FJ analysis should 
be regarded as an important step towards a fuller 
description of fluctuation effects rather than an end in 
itself. Nevertheless the systematic approach which they 
adopted for integrating out degrees of freedom is crucially 
important to subsequent theoretical developments and 
we recall some of the essential features of their formalism. 

FJ seek to derive an intcrfacial Hamiltonian from a 
LGW model and begin their analysis by emphasizing the 
need to carefully define the collective coordinate £(y). 
There is of course a good deal of freedom in the choice 
of definition pTj| (see later) and utility is an important 
criterion. By far the easiest option is to adopt the 
crossing criterion in which ^(y) corresponds to the surface 
of fixed magnetization m x . In their analysis FJ set 
m = but in keeping with subsequent developments 
we keep m x arbitrary. The generalized effective 
Hamiltonian for the surface of fixed magnetization is 
defined as l20ll 



-H FJ [t;m x ] 



Vm e - HLGw[m] 



(1.27) 



c 



where C denotes a constrained functional integral or 
partial trace over configurations m(r) which respect the 
crossing criterion 



i(r = (y,*(y)) 



(1.28) 



Next FJ argue that for a given collective coordinate 
configuration all other fluctuations are small and may be 
ignored. This leads to a saddle point identification 



H FJ [£;m x ] = H LGW [m s (vJ)} 



(1.29) 



where m~ is the magnetization configuration which 
minimises HhGW^A subject to the crossing criterion. 
Fortunately in order to derive effective Hamiltonians 
describing long wavelength fluctuations it is sufficient 
to consider planar constrained profiles m n (z;£ n ) which 
satisfy the Euler-Lagrange equation p4|,p7| 



<(*) = 0'(m*(*)) 
with boundary conditions 
dm n 



dz 



= 0i( TO 7Tl) 



Z = 



lim / \ 

z— >oo ra^yz) = 771,3 



(1.30) 

(1.31) 
(1.32) 



(where m^x — m n (0;£ w )) and the planar constraint 



In this way FJ derive 



(1.33) 



Hit m x ] 



dyl hl^{y)-m x 



V£(y) 




(1.34) 



where the binding potential and position dependent 
stiffness for the surface of fixed magnetization m x are 
determined by m 7T {z;£ v ) as 



W(e(y)-r 

> 

dz 



01 ("Vl) 



1 / dm^ 

2 \~dz~ 



0(m T ) 



(1.35) 



t*=£{y) 

(ignoring some £ independent terms) and 



\ din 



Zap + AE(lfr); 



m 



(1.36) 



Note that the position dependent increment AE — > 
as £ — > 00 for sensible choices m a > m x > mp 
corresponding to surfaces of fixed magnetization which 
unbind from the wall. For these choices the explicit 
expansion of the bind ing potential is very similar to the 
CW expression (1.16) but the AE term is new |pH 



AE(£; 0) = 2nm a Te~ Kl - 2n z m z a VK£e 



(1.37) 



with v = 2— £. Note that by construction, minimization 

of W(£;m x ) recovers the MF position of the surface of 
fixed magnetization m x . 

The (linear) RG analysis of the FJ model is more 
complicated than the CW model because the flows of 
AE and W are coupled. Nevertheless FJ show [0l 22 
that the effective potential (with t = In 6 the usual 
infinitesimal rcscaling factor) 

A 2 

W^{£) = W ( - t \£;0) + ^- Y (l-e- 2t )A^ t \£;0) (1.38) 
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satisfies the same diffusion-type flow equation familiar 
from the linear RG analysis Fisher and Huse made of 
the CW model fl. Due to the negative next-to- leading 
order exponential term in A£(^;0) the critical wetting 
transition is destabilised for sufficiently small lu < lu* and 
is replaced by a weakly first-order phase transition. FJ 
argue that the value of lu* is in the range h < lu* < 1 for 
the Ising model although recent non-linear RG studies of 
the FJ model estimate a significantly larger value lu* m 2 
]38|| . For lu > lu* the critical wetting behaviour is similar 
to the CW model. 

The FJ description of the complete wetting transition 
is the same (at leading order) as the CW model 
because only the first exponential term in (£) is 
required. Thus their model predicts the same value of 
the ads orption critical amplitude as CW theory, that 
is (fL24j). 



E. Three problems of capillary wave theory 

Here we point out three problems faced by the 
standard CW theory of wetting. Two relate directly to 
discrepancies with Ising model simulation studies and are 
associated with fluctuation effects at the upper critical 
dimension d = 3. The third concerns questions of self 
consistency and arises when we compare the CW theory 
with MF studies of correlation function structure in the 
LGW model. We will argue, building on the earlier 
analysis of PB that the final problem provides the key 
to understanding the other two. 

(PI) Ising simulations of critical wetting 

Extensive MC simulations of wetting transitions in 
the d = 3 Ising model appear to show that the crit- 
ical wetting transition is characterised by MF-likc 
critical exponents. Specifically Binder, Landau and 
co-workers [p|JlO|] studied two wetting transitions 
occurring near 0.6Tc and 0.9Tc (corresponding to 
different choices of the surface field) and establish 
the divergence of the surface susceptibility \i = 
^Sr- According to scaling expectations xi should 



dh 

behave as 



Xl =r- 1 X (hT-^) 



(1.39) 



leading to the prediction xi ~ l/iP 2 " 11 along the 
critical isotherm (route (ii) in Fig. [I]). Contrary to 
the CW model the Ising data are rather well fitted 
by the MF result Xl ~ H -1 ^ 2 corresponding to 
i/|| = l. 

Despite widespread discussion in the literature 
only two detailed explanations have been proposed. 
Shortly after the simulation results were reported 
Halpin-Healey and Brezin |3{| used the CW model 
to calculate Ginzburg criteria for various response 
functions. For the surface susceptibility xi they 



show that crossover from MF to non-classical 
behaviour begins when the transverse correlation 
length is close to a value satisfying 



1 



In 1 



£ 2 A 2 



1 



1 + £n A 2 



1 (1.40) 



Typically this yields a value for £|| of a few 
bulk correlation lengths which Halpin-Healey and 
Brezin argued was not quite attained in the 
MC work. However subsequent simulations for 
larger lattices and smaller bulk fields still show 
no substantial deviation from MF behaviour Jl0| . 
Even re-analysis of the data for the susceptibility 
critical amplitude (which is extremely sensitive to 
fluctuation effects) shows only minor deviations 
from classical expectations and is wildly different 
from CW theory Qj. Moreover simulations of a 
discretized CW model shows no discernible MF 
regime and are supportive of the RG predictions 

Hill- 

The second suggestion is the stiffness instability 
mechanism of FJ who argue that the Ising model 
wetting transition is weakly first-order pl| . While 
this may very well be the case (and certainly 
remains a possibility in our theory) on its own it 
does not provide a quantitative explanation of why 
the MC data for xi is MF-like. 

(P2) Ising model simulations of complete wetting 

Recently Binder, Landau and Ferrenberg (BLF) 
JL2 13 1 have studied thin-film Ising models with 
opposite surface fields hi, hp = —hi acting on the 
spins in the planes z — 0, D respectively. The phase 
diagram for this system was originally predicted 
by Parry and Evans on the basis of a MF 
analysis of the appropriate LGW model and shows 
a number of intriguing features. Of particular 
interest is the behaviour of the correlation length 
in the temperature window Tw < T < Tc (and 
h = 0) which probes finite-size effects at the 
complete wetting transition. In the thin-film this 
corresponds to the one-phase regime in which an a[3 
interface sits at the centre of the slab and is subject 
to very large scale fluctuations. The transverse 
correlation length is predicted to be exponentially 
large in the width D reflecting the shallowness of 
the effective binding potential. It is a simple matter 
to construct a CW theory for this system and use a 
linear RG analysis to predict modifications to the 
MF result. This yields @Q 



6|| ~ 



for kD — > oo 



(1.41) 



where is the complete wetting adsorption critical 
amplitude introduced earlier. The prediction of an 
exponentially large correlation length was checked 



G 



and confirmed by BLF for several temperatures in 
the range T\y < T < Tc (it can be extracted 
relatively easily from the the total and mid-point 
susceptibility). However the exponent of the 
exponential term is inconsistent with both MF 
and CW theory! Using their data Parry, Boulter 
and Swain (PBS) |Q have extracted the effective 
value of the critical amplitude 9 as a function of 
temperature. A plot of the simulation results is 
shown in Fig. || and shows a substantial increment 
to the expected value 9cw ~ 1-4 (recall 9mf ~ 
1) for temperatures deep in the complete wetting 
regime (T > T w » 0.97b). 
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FIG. 3. BLF data for vs. T /Tc showing the anomalous 
enhancement to the CW result 6 — 1 + u/2 m 1.4 deep 
in the complete wetting regime. Apart from the point at 
T /Tc = 0.916 the error is within the symbol. The dotted line 
is a cubic fit to the data while the asterix is the extrapolated 
value of 8 at Tw ~ 0.9Tc- The decrease in the value of 9 
as T —* T w is indicative of novel crossover behaviour near 
critical wetting associated with coupling effects. 

Specifically at high temperatures the simulation 
results are consistent with 9 w 1.7-1.8. The 
abnormally large values of 9 are inexplicable 
using the CW and FJ models and raise further 
doubts whether the standard picture of interfacial 
fluctuation effects captures all the essential physics 
at the upper critical dimension. 

(P3) Correlation function structure 

While the importance of the two problems de- 
scribed above can be recognised immediately it is 
not at all clear how to modify existing effective 
Hamiltonian theory to overcome them (assuming 
that this is where the fault lies). However there 
is a third long-standing puzzle associated with 
wetting theory, which while less dramatic than 
the questions concerning non-universality, focuses 



attention on a specific defect of any model based 
on a single collective coordinate £(y)- 

MF studies of the order parameter correlation 
function G(ri, r 2 ) for a number of different density 
functional models of wetting reveal intriguing 
features for positions close to the wall which 
cannot be described using the CW model (or for 
that matter the amended model proposed by FJ) 
|44[|45|. It should be emphasised that the MF 
studies are consistent with some curious exact sum 
rule requirements which link correlation function 
and thermodynamic singularities and which also 
defy explanation using standard arguments |45| , [46f . 
We expect therefore that much of the correlation 
function structure described below is also present 
in the full three dimensional LGW model beyond 
MF approximation. 

First, consider the complete wetting transition in 
the LGW model. The MF correla tion function 
satisfies the differential equation (1.21) whose 
solution has a simple Lorentzian form for z\ , z-i near 
the a[3 interface. However for particle positions at 
(or close to) the wall the solution is non-Lorentzian 



G(0,0;<?) 



~ 'i 



(c + K)m'i + q 2 



(1.42) 



where fh! x is the gradient of the magnetization at 
the wall (which is of order r and may be regarded 
a constant). As emphasised by PB the wave- 
vector dependence shows crossover from coherent 
to intrinsic behaviour depending on the scaling 
variable x = q£u . By this we mean that for 
finite £|| and q — ► the effective coefficient of 
q 2 in the denominator is the full surface tension 
(or more properly, stiffness) of the wall-/? interface 
and indicates the coherence of asymptotically long 
wavelength fluctuations in the wetting film. This 
manifests itself in the second moment, defined in 



the expansion G(zi,Z2;q) 
such that 



G 2 (0,0) tx a w /3 - <£i (mi) 



(1.43) 



This is no artifact of the MF approximation since 
an equivalent statement is known to be precisely 
true for the phenomenon of complete drying by a 
vapour phase (/3) at a hard wall-liquid (a) interface. 
The exact sum-rule reads [M 



— G2(0, 0) — a wa + a aj 3 + / s j ng 
and is valid for arbitrary fluid-fluid forces. 



(1.44) 
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On the other hand in the scaling limit of finite 
q and £|| — > oo (i.e. h — > _ ) the correlation 
function reduces to the appropriate expression for 
fluctuations typical of the wall-Q interface. For 
this case the effective coefficient of q 2 in the 
denominator is the local, intrinsic stiffness H wa = 
o wa - 01 (mi). 

The correlation function with one particle at the 
wall and one at the a[3 interface also shows 
interesting features. For example in their MF study 
of drying at a hard-wall within the Sullivan density 
functional model Parry and Evans derive [H 



G 2 (0,l) 
m'{£) 



" k m 



(1.45) 



where here m'(£) denotes the gradient of the 
equilibrium number density at the vapour-liquid 
(a/3) interface z — £. The second and third terms 
on the r.h.s. correspond to logarithmic next-to- 
leading order and coherent properties, respectively. 
The same behaviour emerges in the LGW model at 
MF level although, unfortunately, there is no exact 
sum-rule to compare with. 

A cause for concern for theorists is that hardly 
any of these features are describable using one-field 
Hamiltonians such as the CW and FJ models. In 
fact one of the virtues of the FJ approach is that one 
may unambiguously calculate a correlation function 
G FJ (zi,z 2 ;q) (for a given m x ) at MF level 
(and beyond) since the theory allows a one-to-one 
connection between order parameter and collective 
coordinate configurations through the field £(y). In 
this way it is straightforward to show that the FJ 
model with m x = (and by implication the CW 
model) cannot recover the non-Lorentzian, next- 
to-leading order and coherent behaviour manifest 
in the MF G(zi, z 2 ; q) at complete wetting [2j|. 

Similar remarks apply to the critical wetting tran- 
sition. Most importantly the correlation function 
the wall — which can be 
noting that rh'i ~ r and 
The result is conveniently 



is non-Lore ntzia n 
read from ( 1.42 ) 
v wa ~ 0i (mi) ~ 
written 

1 



at 
on 

_2 



1 



const. 



G w P(0, 0;q) G wt *(0,0;q) 1 



(1.46) 



and shows the necessary crossover from singular 
to intrinsic behaviour dependent on the scaling 
variable x = . This expression is consistent with 
a number of sum-rule results which establish precise 
connection between the moments of G(0, 0; q) and 
thermodynamic singularities |^,^|. The most 
important of these are 



G o (0,0)~7 
G 2 (0,0)~a Q/3 T 



and recall that a s = 
addition we note that 



dzG o (0,z) = 



in MF theory. In 



dm i 
~dh 



(1.48) 



-2(1+/3 S 



(1.47) 



which shows the role played by correlations from 
the interface to the wall. 

A more precise MF expression for G w ^(0,0;q) 
at critical wetting will be derived later using the 
generalized stiffness matrix formalism. The point 
that we wish to emphasise here is that the CW 
and FJ models do not satisfy the exact sum-rules 
in a manner that is also consistent with the non- 
Lorentzian form of G(0, 0; q) seen in the MF LGW 
calculations. The reasons for this will become 
clearer in the next few sections. 



F. Outline 

In this paper we further develop the idea forwarded 
by PB that the problems of wetting theory reflect the 
failure of effective Hamiltonian models based on a single 
collective coordinate to account for the coupling of 
interfacial fluctuations. Thus the physical mechanism 
underlying our mathematical modelling is the following 
coupling hypothesis 

The correlation function near the wall and 
some of the non-universal physics associated 
with wetting at the upper critical dimension 
d = 3 is sensitive to the coupling of order 
parameter fluctuations near the unbinding 
interface and wall. 



PB concentrated on the complete wetting transition and 
constructed a two-field model H\£\,&2\ which in Gaussian 
appro ximation (generalising the calculations leading to 
( 1.20| )) identically recovered the MF correlation functions 
near the interface and wall. Here the generalized 
collective coordinates i\ and £2 represent surfaces of 
fixed magnetization which remain bound and unbind 
from the wall respectively. Subsequent RG ]l6| , p7| , ^9[ 
analysis of the two-field model for fluctuation effects in 
d = 3 showed that the effective value of the wetting 
parameter determining the critical amplitude 9 was 
renormalized compared to the CW expression. Thus 
addressing (P3) for complete wetting appears to resolve 
the second problem (P2), at least semi-quantitatively. 
Unfortunately the PB model is not well suited to 
studying critical wetting and only hints at a possible 
influence of coupling on this transition Q . 

The task of the present work is to generalise the PB 
approach to obtain a quantitative theory of coupling 
effects at complete and critical wetting. To do this 
requires a more thorough analysis of the connection 
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between microscopic models (such as the LGW Hamil- 
tonian) and collective coordinate theories. As we shall 
show it is first necessary to broaden the generalized 
effective Hamiltonian approach to allow for different 
types of collective coordinates describing order parameter 
fluctuations near the wall. Having done this we argue 
that it is possible to chose an optimal definition using 
a novel variational principle such that the sum over 
all possible configurations appearing in the effective 
theory most accurately mimics the partition function 
for the microscopic model. The optimal Hamiltonian, 
written H[s,i], is then studied to elucidate the influence 
of coupling on the surface phase diagram and critical 
properties. 

Our paper is arranged as follows: we begin by 
reviewing the PB theory and its successes in addressing 
(P2) and (P3) for complete wetting. The breakdown of 
the PB description near critical wetting is then discussed 
in addition to some pointers for possible coupling effects 
at this transition. In Sec. HI we introduce further 
examples of coupled Hamiltonians which use different 
types of collective coordinate at the wall. We provide a 
proof of a generalized correlation function reconstruction 
scheme (CFRS) which shows that the physical correlation 
function G(zi,Z2;q) may be recovered (at MF level) 
using arbitrary choices of the collective coordinates, 
provided they satisfy certain conditions of locality. The 
invariance of G{z\ , z 2 ; q) is crucial to our further analysis. 

In Sec. |fv| we introduce a particular class or set Ti 
of coupled models. This smoothly interpolates between 
two cases which describe order parameter fluctuations 
near the wall as local translations and enhancements 
of the magnetization respectively. The set 7i contains 
more general types of collective coordinate (which we 
shall refer to as proper coordinates) which have both 
interfacial- and spin-like components. All of the models 
in H are candidates for describing coupling effects beyond 
MF theory. We forward a novel optimization scheme 
leading to a unique choice of proper collective coordinate 
model in Ti which best describes fluctuations at the wall 
coupled to the unbinding interface. 

In Sec. IV B we consider fluctuation effects and show 



that the extent of the coupling is controlled by an angle- 
like variable S* whose value depends on the proximity 
of the critical wetting phase boundary. In particular 
for T 3> T\v deep in the complete wetting regime the 
angle 6* ~ § and the optimal theory reduces to the 
PB model. As the temperature is reduced, so that 
the thermodynamic path (iii) is closer and closer to 
TV, the angle 5* rotates to zero signifying a qualitative 
change in the role of coupling. This mechanism yields 
a precise expression for the temperature dependence of 
the renormalized wetting parameter determining 9 in the 
proximity of T\y- For critical wetting the coupling has 
the effect of dramatically decreasing the extent of the true 
asymptotic critical regime for the local susceptibility Xi 
consistent with the Ising model MC observations of MF- 
like behaviour (problem (PI)). The calculations for both 



critical and complete wetting highlight the role played 
by a second wetting parameter (written fl) associated 
with fluctuations near the wall which has no counterpart 
in CW and FJ theories. We discuss the temperature 
dependence of the new wetting parameter and argue 
that it approaches a universal value Qc close to unity 
as T^Tp. 

The justification for the coupling hypothesis is largely 
a posteriori and relies on the success/failure of the 
theory to explain problems (P1)-(P3). Nevertheless we 
conclude our article with some remarks suggesting that 
the coupling to order parameter fluctuations at the wall 
can only effect thermodynamic singularities at the upper 
critical dimension for systems with short-range forces. 



II. STIFFNESS MATRIX FORMALISM FOR 
COMPLETE WETTING 

We begin our discussion of coupled fluctuations by 
reviewing the formalism and results of the PB theory of 
complete wetting which will be generalized in subsequent 
sections. To this end we first summarise the main results 
establishing the connection between FJ theory and MF 
correlation functions (further details may be found in 



A. FJ theory and correlation functions 

The systematic structure of the FJ theory allows 
precise connection to be made with the MF correlation 
function of the L GW theory satisfying the Ornstein- 
Zernike equation (1.21). Following Parry p5| , consider 
the continuous set of FJ Hamiltonians {Hpj\£; m x ]} by 
allowing for all possible values of m x belonging to the 



range of magnetizations rh\ > 



> mp seen in the 



MF profile m(z). Ignoring fluctuations, the equilibrium 
position £ = z of the surface of fixed magnetization 



clearly satisfies m 



x 



fh(z). Now consider the 



corresponding structure factor 



5(q; z) = J dy 12 e^((£( yi ) - z){i(y 2 ) - z)) (2.1) 

In Gaussian approximation, valid for small fluctuations, 
this is trivially calculated 



S(q;z) 



W"(z;m(zj\ +q 2 z(z;m(z) S j 



(2.2) 



where 



W"(z;m{z)) = -^W(£;m x ) (2.3) 



for I — z and m x = fh(z). Note that in writing (2.2) 
we have ignored terms of order q A that would arise if 







we included an expression for the rigidity in the FJ 
model. However these are not important for interfacial 
phenomena and will not be discussed further (see Ej| 
for explicit calculations). It transpires that there is a 
remarkably elegant relation between the set {S(q;z)}, 
parameterized by z, and the MF correlation function 
G(z, z; q). To see this note that within FJ theory a small 
local translation (say) in the position of the surface of 
fixed magnetization m x generates a precise change in 
the magnetization 



5m(y) 



(z;£^S£(y) 



(2.4) 



at position z. Therefore a Gaussian approximation 
for each of the FJ Hamiltonians in the set may be 
used to calculate an expression for a magnetization 
pair correlation function at any two points. Let us 
denote the result for the Hamiltonian H[£;m x ] by 
G FJ {z\, z 2 ; q; m x ). Then we can identify 



G J (zi,z 2 ;q;m 
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{z 1 ;e^)-^-(z 2 ;£ 7r )S(q;z) 



(2.5) 

where the values of m x and z on the l.h.s. and r.h.s. are 
related by m x — fh(z) and the partial derivatives of the 
planar constrained profile are evaluated at equilibrium 
£ Tr = £ = z. The question now remains how does this 
set of correlation functions relate to the actual s olutio n 
G(z\, z 2 ; q) of the MF Ornstein-Zernike equation (p.. 21 )? 
Explicit calculation shows that each element Hpj[£; m x ] 
generates the correct MF correlation function only for 
positions z\ — z 2 — z. Specifically we can identify 



G(z, z; q) = G FJ (z, z; q; fh{z)) 



which becomes 



G(z, z; q) = rh'(z) 2 S(q; z) 



(2.6) 



(2.7) 



on using 



d£„ [t ' 



-m'{£) 



(2. 



(obta ined by simple differentiation of the crossing crite- 
rion ( 1.2§| )). For example, explicit evaluation of (2.2) in 
FJ theory yields 

G (z,z) =m'(z) 2 S(0;z) 
= m (z) 



fh'i (cfh'i — m") Jo fh'(z') 2 



dz' 



(2.9) 



which is the analytic solution to (1-21). These remarks 
make it clear why the FJ model with m x = (and 
hence the CW model) fail to describe the wall correlation 
function G(0, 0; q) since the choice of m x is entirely 
inappropriate. With m x = the FJ model accurately 
describes fluctuations near the a(3 interface and can 
exactly recover the MF G(z, z; q) only at the position 
z where m(z) = 0. 



B. Two-field theory and correlation functions 



PB next consider the properties of two-field models 



H[£ u t 2 y, 
and m x 



,A"1 



X 



for surfaces of fixed magnetization m\ 
described by a pair of collective coordinates 



£\ and £ 2 ||15j. Importantly for the complete wetting 
transition, on which they concentrate, the MF and FJ 
profile has a significant lip at the wall (see Fig. ||) so that 



it is possible to chose fh\ 
so that 



> 



> mp and m x — (say) 



(i) the upper surface unbinds 

(ii) the lower surface remains close to the wall 
as ft. 0. 

magnetization 




m It (z;^ 2 ) 



FIG. 4. Schematic portrayal of the single (dashed line) and 
doubly constrained planar profiles near a complete wetting 
transition. Discontinuities in the gradients exist at z = li 
and z — £2 but these have not been explicitly shown. 



Consequently the MF positions £ fl 



with ji = 1, 2 



of these generalized surfaces represent planes that are 
close to the wall and unbinding interface. The two-field 
Hamiltonian is constructed using a natural generalization 
of the FJ method. The fundamental saddle point 
identification is 



H[t x 



= H 



LOW 



[r4 2) (r;4,4)j 



(2.10) 



where we have dropped the explicit m x dependence. 

(2) 

Here m~ represents the magnetization distribution 
which minimises the LGW Hamiltonian subject to the 
double crossing criterion 



x 



withr=(y,^) (2.11) 



As mentioned earlier calculations are eased when we 

observe that only the properties of planar constrained 
(2) 

profiles m\ '(z;£i,£ 2 ) need be considered to derive 
models with gradient terms. 
In this way PB find jl5| 

(2-12) 
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where the binding potential is 



It 



dz 




<t>( m n 



(2) 



(2.13) 



ignoring some constant terms independent of £ 2 (which 
are subtracted off). The stiffness matrix elements are 
given by 



£^(4,4) 



(2) 



dz 



d£ 



/ITT 



BP 



(2.14) 



i(y) 



The constrained profile satisfies the same Euler- 



Lagrange (1.30) considered earlier but must now be 

solved subject to the double planar crossing criterion. 

(2) 

Once (z;£i,£ 2 ) is found then the explicit position 
dependence of W 2 {£\,£ 2 ) and £^(-^1,^2) may be calcu- 
lated. Conveniently the binding potential decomposes 
into a shielded form 



W 2 {£ u £ 2 ) = U(£ 1 ) + W(£ 2 -£ 1 ) 



(2.15) 



where W{£) is similar to the CW result (1.16). The 
first term simply serves to bind the lower surface to 
the wall. We will always assume that the fluctuations 



-1/2 



ee [{8£\)_ 
complete wetting 



are small so that we may write for 



W 2 {h,£ 2 ) = -£\ 



h£ 2 \ + 2Km a Te 



-Ki21 



(2.16) 



with r ~ r 2 explicitly determined and 



21 = h- h- 

Minimization of W 2 (t\,£ 2 ) gives the MF positions z\, z 2 
of the collective coordinates. The position dependence of 
the stiffness coefficients can also be determined 



En 



<J a p + 0{e~^) 



£12 ~ n m a rf 2 ie" 



(2.17) 



and notice that a wa — 4>i{fhi) ~ r 2 . In these expressions 
fh\ is to be interpreted as the magnetization for 
the wall-a interface. In fact the position dependence 
of the stiffness coefficients plays no part in determining 
leading order thermodynamic effects although the cross 
term £12 is important as regards the coherent and next- 
to-leading order correlation function behaviour discussed 
earlier 

Due to the inclusion of two collective coordinates local 
to the a/3 interface and wall, the coupled model can 
fully explain the non-Lorentzian structure of G(0,0;q) 
described in (P3). In fact precise connection with MF 
correlations can be made using the following CFRS p8) . 



Define the structure factor matrix S(q) as the 2x2 matrix 
with elements 



S^q) = J dy 12 e l - 



! (^(yi)^(y 2 )) (2.18) 



These can easily be calculated using the relation 



S^q) = ( S 1 d h 



r) 2 ft 2 
u 21 °22 



W{t u £ 2 )+q 2 V 



(2.19) 



where d^ v 



d- 



and is evaluated at equilibrium. 



From the matrix elements one can calculate the MF 
expressions for three possible correlation functions by 



(2.20) 



where z^Zv e 



{zi,z 2 }, i.e. — m(z M ). This 

identically recovers the known solution to the Ornstein- 
Zernike equation (1.21). 

For example, at the wall PB derive 



™'2 



G(0,0;q) = 



m 



En 



(2.21) 



•0(g*) 



which is of the required form fll.42j ) . PB proceed to show 
that the elements of the stiffness matrix obey the relation 

US 



^ £ M „(0, 22) = CTwa - <f>i(mi) 



(2.22) 



which ensures that the sum-rule (1.43) describing the 
coherence of long wavelength fluctuations is exactly 
satisfied. The origin of the / s i ng contribution to G2(0,0) 
is elegantly explained by this approach through the 
position dependence of the off-diagonal stiffness matrix 
term since 



2£i 2 (0, Z 2 ) ~ /sing 



(2.23) 



and recall that the arguments refer to the mean positions 
of l\ and l 2 . Notice that the position dependence of 
£12 ( see ( [2-17 )) is longer ranged than £ n and £ 2 2 and 
is precisely of the form required to yield the singularity 
/sing ~ 7z In 1 7i I . Following PB we refer to ( 2. 22] ) and 
( |2.23 ) as stiffness matrix free energy relations. 



C. Renormalization of the wetting parameter 

The most important consequence of coupling is inde- 
pendent of the position dependence of the stiffness matrix 
elements and it suffices to write the Hamiltonian [fl7| 



H[£ 1 ,£ 2 ] = H [£ U £ 2 ]+H 1 [£ 1 ,£ 2 



(2.24) 



where 



11 



and 



1, 



1 



dy i -Zu(WiY + -£ 22 (W 2 ) 2 + -rl{ 



Hi[£i,e 2 ] = / dyW{£ 2 ,£ x ) 



(2.25) 



(2.26) 



In this approximation we may identify E 22 = £«/3 
and En = o wa — 0i(mi). Note that \J^^- = £,wa 
corresponds to the finite transverse (second-moment) 
correlation length for surface correlations at the wall- 
a interface. This can be reliably calculated in MF 
approximation away from Tq yielding £ wa — [k(c + 
k)] -1 / 2 . The momentum cut-off A 2 for the upper field 
is the same as that for the CW and FJ models (denoted 
A). Similarly the momentum cut-off Ai for the lower 
field must be chosen appropriately such that the r.m.s. 
fluctuations are small (and typically of order a bulk 
correlation length) which requires Ai -C V^n f^jl - 

As the fluctuations of £\ are adequately described in 
Gaussian approximation it is possible to account for their 
renormalization exactly. The interaction H\ [£\ , £ 2 ] is 
then treated in linear approximation analogous to the 
RG theory of Fisher and Huse §§ and FJ @. First the 
fast modes (appearing in the Fourier representations of 
the fields) with wave-vectors q satisfying A M /& < q < A M 
are integrated out and the Hamiltonian rescaled leaving 
H [£i,£ 2 ] invariant. The flow equation for W w (£ 1 ,£ 2 ) 
reads 



2 - 







d 2 



- 0J2 



W {t) {£i,£ 2 ) = (2.27) 



dt 1 ~ x d£\ 

with t = In b the infinitesimal rescaling factor and 



K 



U>2 



1 



47tEh \ SuAf + r 
n 2 



4ttS 



a/3 



(2.28) 
(2.29) 



The diffusion equation (2.27) can be easily solved 



■21 



4irt^/ujiuj 2 



d£[d£' 2 W^ (£[,£',] 



<MM- ( ^? 2 - ( ^:? 2 ) (2-30) 



Au> 9 t 



where (£i,£ 2 ) is the bare interaction contribution to 
the binding potential. In the PB model the fields only 
interact through a relative term which is a function of 
£21 = £2 — £\ so that ( 2.30| ) simplifies to 



■y/47r(wi + u 2 )t 



d£'W {0 \£')e~^i+^n 
(2.31) 



which needs to be studied at the matching point t* || 
where the curvature ^ {{£21)) = ^ a pK 2 - At this 
point we may identify e as k£u . These remarks establish 
that the binding potential renormalizes exactly as in the 
CW theory but the divergence of {£) is determined by a 
renormalized wetting parameter 



UJ\ + UJ 2 



4ttEii 



1 - 



EuAf 



(2.32) 



The value of ui lies between two extremes corresponding 
to r = 00 and r — representing cases where the £\ 
fluctuations are suppressed and free respectively. 

The implication of the PB model is that the adsorption 
critical amplitude 9 is increased due to coupling effects 



1 



4 7 r£ii[l + (Ai£ tUQ )- 



(2.33) 



which should be compared to the CW result (1.24) 
(assuming that u,Q < 2). PB estimate Q that the 
increment A9 to the CW result is about 0.3 for Ising- 
like systems deep in the complete wetting regime yielding 
a value for 9 1.7 close to the BLF simulation result. 



Importantly the linear RG expression (2.33) is supported 
by non-linear RG analysis [jl9| and also by MC simulation 
of a discretized version of the model [p0| . 



D. Interpretation 

Before we discuss the unsatisfactory features of the 
PB model we make some new remarks conc ernin g the 
i nterp retation of the fundamental results (2.32) and 
33). The first of these is that the renormalization of 
can be understood heuristically 



(I 



the critical amplitude 
using a simple generalization of the elegant scaling 
analysis of Lipowsky and Fisher [Q which takes into 
account the asymptotic coherence of long wavelength 
fluctuations at complete wetting. 

To see this we initially recover the standard result 

( 



1.24 ) using scaling arguments based on the CW model 
). Neglecting fluctuation effects the singular 
contribution to the free energy per unit area A of the 
wall is clearly 



f s f/ g = W(£) 

hln\h\~ 



0(h) 



(2.34) 



corresponding to MF theory. Now following Lipowsky 
and Fisher we suppose that the fluctuations in the 
position of the a(3 interface lead to an extra bending 
energy so that 



f ■ 

./ SI 



CW 
sing 



1, 



+ f MF 
1 J sing 



(2.35) 
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Using the known interfacial roughness (l.£) and the 
diverg ence of £|| (which is unchanged from the MF result 
(1.14)) we are led to 



/, 



cw 

sins 



,„ w. hhx \h\~ 



(2.36) 



which is the explicit RG result (for uj < 2) consistent 
with flL2§) (recall that {£) = ^ E ). 

To include coupling effects we first assume (as is 
reasonable) that the estimate of the bending energy 
contribution from the large scale capillary-wave-like 
fluctuations of the ^ 2 (y) field is unchanged i.e. the 
roughness relation between £± and £ii of the a/3 interface 
is the same as the CW model (see below). Further, we 
suppose that the only effect of the small fluctuations in 
the £i{y) field is to increase the effective area of the wall 
as seen by asymptotically long wavelength (coherent) 
fluctuations in the a(3 interface. The relative increase in 
area can be estimated (analogous to the bending energy 
contribution) as 



where 



K>) ; 



k R T K 2 



47rEn 



In 1 + (Ai& 



(2.37) 



(2.38) 



determines the r.m.s. fluctuation of the lower surface. 
Now due to the asymptotic coherence of long wavelength 
fluctuations the direct MF contribution to the singular 
free energy must be modified due to the increase in the 
effective area of the wall. Thus we must add a coupling 
contribution 



W(£) (2.39) 
which with the usual MF bending energy term gives 

/£«fl + S + ^(AxW a> )^^ (2-40) 



and correctly identifies the renormalized expansion for 9 
for small K\. 

A nice feature of this argument is that it directly 
relates the renormalization of to the coherence of 
asymptotic long wavelength fluctuations signifying the 
importance of similar coupling effects for problems (PI) 
and (P3). 

Implicit in the above is the assumption that the 
roughness of the a [3 interface is not changed by coupling 
effects. This takes us to our second remark that 
statements concerning the renormalization of the wetting 
parameter only refer to the critical amplitude 6 and 
not the value of uj determining the divergence of the 
roughness £j_ • In fact this can be easily deduced from the 
renormalized CFRS which identifies the structure factors 



Suv(q) (defined in ( 2.1 8| ) ) beyond MF approximation in 
d = 3 g|. 

The renormalized CFRS is very similar to that 
described earlier because the RG matching procedure 
reduces the problem to a Gaussian calculation just as 
in MF theory. Thus we need only replace the binding 
potential and stiffness matrix by their renormalized 
values and allow for the spatial rescaling factor. We first 
calculate the renormalized (inverse) matrix 



( S ( t *)(q) 



r- 








U 2 i 





Wp{h,£ 2 ) + q z ^e 



21* 



(2.41) 



where 

=e 2t ^r£l + W^(t2i) 



„2J 



r£\ + h£%\ + 2m a KTe' 



(2.42) 



and we have assumed that uj < 2. Ignoring the position 
dependent elements, the stiffness matrix is 



En 
Iw 



(2.43) 



and is not renormalized (analogous to the stiffness 
coefficient in CW theory B). The required structure 
factors then follow as 



(2.44) 



which are easily calculated and yield final results that are 
basically unchanged from the MF calculations described 
earlier. In particular the 522(9) element determines the 
mean square value of the Fourier amplitude 8 £2 (q) (see 
PI) 



S2 2 (q) = (|«2(q)| : 



JSt* 



-,21 



*E a/3 g 2 



W(** )" +e 2t * (r+g 3 Si 1 ) 



•21* 



WW + E a/ 3 9 2 e 2t * 
1 

E a/3 (^ 2 + q 2 ) 



(2.45) 



using -* = S Q /3K 2 and e* = Performing the 

Fourier inversion in the standard fashion yields 

(< ± ) 2 ^ln(l + (A£||) 2 ) (2.46) 

as quoted in the introduction for the CW model. Thus 
the effective value of the wetting parameter is unaltered 
as regards the interfacial roughness relation in contrast 
to the calculation involving the critical amplitude 8. 
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Combining ( 2.46| ) with the two-field prediction for the 
transverse correlation lengt h in the thin-film geometry 
with opposite surface fields (1.41) yields 



for D — ► co 



(2.47) 



in the soft mode regime. This shows that the divergence 
of the roughness as D — ► oo is controlled by the bare 
and renormalized wetting parameters (see ( 2.33| )). Data 
for the mid-point width (denoted asra = \/2Tr£,± 

in [pl|) have been extracted from the BLF simulation 
data |51| at ^- ~ 0.9554 and compare favourably with 
the asymptotic prediction above for moder ately wide thin 
films D < 40. Note that the prediction ( 2.47 ) correctly 
identifies the pre-factors for the square root divergence 
of £_l (and hence w) misquoted in |5l[] although the 
numerical estimates for w are very similar to their Fig. 
3(b). 

Also note that the present RG result for Su(q) gives 
an identical non-Lorcntzian expression as that found in 
the MF calculation ( 1.42] ) strongly suggesting that the 
exact G(0, 0; q) has such a structure. 



E. Problems with the description of critical wetting 

1. Crossover from complete to critical wetting 

The PB model is certainly not capable of yielding a 
fully quantitative theory for the temperature dependence 
of the renormalized critical amplitude 9 which could be 
compared in detail with the simulation results. This is 
due to the uncertainty in the value of the momentum 
cut-off Ai for the lower field which appears in the 
RG prediction. This is a serious defect of the model 
which forces us to rethink the nature of order parameter 
fluctuations near the wall in the vicinity of T\y- To 
see this we follow PBS and consider complete wetting 
transitions occurring along paths (iii) in Fig. [I] closer 
and closer to the wetting temperature TV. Recalling 
that the stiffness coefficient for the lower surface vanishes 
like En ~ t 2 we are forced to conclude that the cut- 
off Ai must also vanish in th is limit otherwise the value 
of Q is singular (see ( 2.33 )). In fact such behaviour 
can be anticipated since as mentioned earlier we need 
to impose the inequality Ai <C %/En in order that the 
r.m.s. amplitude of the Gaussian fluctuations ^i(y) at 
the wall remains bounded. This pathology is intimately 
linked to the behaviour of the MF and FJ profiles near 
the wall which flatten as the critical wetting transition 
temperature is approached, signalling the breakdown 
of the crossing criterion method used by PB to model 
small order parameter fluctuations. In fact this feature 
is also seen in the equilibrium magnetization profiles 
extracted from the Ising model simulation studies — see 
for example Fig. 4 of |13|. PBS point out that this 



has important consequences for complete wetting close 
to Tw since it implies that the increment to u> due to the 
coupling of fluctuations must vanish 



so that 



lim ( 



i 



(2.48) 



(2.49) 



This prediction is in very good agreement with the BLF 
simulation data (see Fig. 0) and extrapolation to T\y ~ 
0.9T C yields 9 w 1.4 and hence w(T w ) w 0.8 which 
is close to the series expansion estimate of Fisher and 
Wen at th is tem perature |j36| . As conjectured by PBS 
the result ( 2.48 ) is indicative of some kind of decoupling 
phenomenon between order parameter fluctuations near 
the a[3 interface and wall occurring close to Tw ■ However 
the PB model can only describe this in the crudest 
possible (and unsatisfactory) way by allowing the cut-off 
Ai to vanish as T — > T^. 



2. Towards a theory of coupling at critical wetting 

Despite the unwelcome behaviour of the cut-off Ai 
during the approach to T\y we note that the two-field 
theory does hint at a possible explanation of problem 
(PI). As pointed out by Parry and Boulter |4g] if we 
simply assume that the model, with a suitably chosen 
cut-off Ai <C t 2 , is appropriate for describing coupled 
fluctuations along route (i) in Fig. [l] then it immediately 
follows from the RG transformations described earlier 
that the structure factors Sn(0) and Si2(0) associated 
with the lower field retain MF-like singularities even 
whe n S22 (0) has a non-classical divergence with v\\ given 
by ( 1.22 ). While we should be rightly suspicious of the 
precise status of this prediction one might tentatively 
suggest that coupling somehow reduces the manifestation 
of critical singularities for local response functions at the 
wall. As we shall see in the next sections this conjecture 
is supported by a more quantitative theory of coupling 
effects at wetting transitions. 



F. Open questions 

While the semi-quantitative successes of the coupled 
model introduced by PB are encouraging there are a 
number of unpleasant features associated with the theory 
that warrant further investigation. All of these concern 
the pathological behaviour of the cut-off Ai associated 
with the approach to the wetting temperature which 
clearly limit the application of the model. This is most 
unwelcome from a fundamental point of view since we 
should take a harsh opinion of any field theory (at least 
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in condensed matter physics) in which the cut-off has 
singular behaviour. 

Our central goal in this paper is to develop a coupled 
theory of wetting transitions which involves a unified 
description of order parameter fluctuations at the wall. 
It should be equally applicable to complete and critical 
wetting and avoid the need to impose any special 
behaviour of the cut-off (s). The theory should also be 
able to explain the apparent success and simplicity of 
the PB model deep in the complete wetting regime and 
provide a more sensible mechanism for t he de coupling of 
fluctuations associated with the result (2.48) (assuming 
that this survives). 



III. GENERALIZED COUPLED MODELS AND 
THE CFRS 

To develop a fully quantitative model of coupling 
effects at wetting transitions it is necessary to reassess 
the role played by collective coordinates in effective 
Hamiltonian theory. The original motivation of PB was 
to construct an effective Hamiltonian which in Gaussian 
approximation could reproduce the Lorentzian and non- 
Lorentzian structure of the MF correlation function near 
the af3 interface and wall respectively. As PB showed 
MF correlations could be precisely recovered using the 
CFRS described in the last section. We shall now that 
show that this is also possible for different choices of 
the collective coordinates. We begin with an illustrative 
example before considering the general case. 



A. Interfacial- and spin-like collective coordinates 

In Sec. |n| we reviewed how the two-field theory 
grew out of observations concerning the connection 
between MF correlations and the set of FJ Hamilto- 
nians {H[£;m x ]}. Heuristically, the crossing criterion 
naturally defines an interfacial-like collective coordinate 
£(y) which generates magnetization configurations by 
translating the position of a surface of fixed magneti- 
zation. However it is also possible to define spin-like 
collective coordinates which model local enhancements 
of the magnetization in some plane. We begin by 
considering one-field models parallelling the analysis of 



Sec. [I A 



Let us focus on the behaviour of the d—1 dimensional 
plane of spins a distance z from the wall and constrain 
them to satisfy 



(3.1) 



We may define an effective Hamiltonian for this spin 
plane by tracing over all magnetization configurations 
that leave {o(y)} invariant 



H[a; z x ] = - In 



Dme" fltGwIm| 



(3.2) 



where C denotes the spin constraint (|3.1| ). Now let us 
assume, rightly or wrongly, that all fluctuations apart 
from those of the spin plane arc small such that we can 
use a FJ-like saddle point identification 



H[a;z ] = min H LGW [m] 
= H LGW [m-s] 



(3.3) 



where tyi*?(y', o ) is the profi le t hat minimises Hlgw^kA 
subject to the constraint (3.1). Whether or not this 
assumption is justified will be considered later. For our 
present purpose of reconstructing MF correlations it is of 
no consequence as we now show. 

The constrained magnetization satisfies the usual 
Euler-Lagrange equation 



V 2 m= + </>'(m s ) = 



(3.4) 



with boundary conditions 



dz 



Mm I \ 

z—>oa fflB^J = 171/3 



(3.5) 



and spin criterion (3.1). 

In the same manner as FJ theory however it suffices 
to consider planar constrained profiles m„(z; ov) from 
which ma(r; o) may be constructed perturbatively using 



/ \ Qvn 

m H (r; a) = m ff \z; a(y)J + -^-Sa(y) 



(3.6) 



with 5a(y) = o(y) — ov. In this way it is straightforward 
to construct the effective spin plane Hamiltonian 

H[a;z x } = J dy^(a;z x ) + W(a;z x )^ (3.7) 

(ignoring terms of 0(V 2 o)) where 



W (a; z x ) = i<j>i{mAQ; <**)) 



dz 



1 / dm^ 

2 \~dz~ 



(r»v(z; ov)) 



0v=<r(y) 



and 



(3.8) 



(3.9) 



<r*=<r(y) 



By construction minimisation of W(a;z ) recovers the 
MF value of the magnetization at position z — z x . Thus 
we can identify 
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with 



rh(z 



(3.11) 



So far our analysis parallels that of FJ but with a 
different type of collective coordinate. Now we follow 
Parry |2f| and use the set of spin-plane Hamiltonians 
(by allowing all possible values of z x ) to derive a set of 
correlation functions G^'(zi, z 2 ; q; z x ). Consider first 
the zeroth moments G^ (zi, z 2 ; z x ). Each Hamiltonian 
can be used to calculate a correlation function since there 
is a one-to-one relation between fluctuations in a and 
the magnetization at any point v\. For homogeneous 
enhancements we note that 



5m(ri) 



<9ov 



(zi;ov)<5cr(yi 



and in particular at the plane z 



itself 



dc 



(z x ;a) = l 



(3.12) 



(3.13) 



Assuming that fluctuation are small and can be ade- 
quately described by a Gaussian approximation we arrive 
at the prediction 



G^ ) (z 1 ,z 2 ;z x ) 



-Q^-{z 1 ;a 7T )- 5 ^-{z 2 ;a 7T ) 



W"(a;z x ) 



(3.14) 



for the (zeroth) moment correlation function according 
to the model (|3.7[). Here 



W"{a;z x ) = 



d 2 W 
~d^ 



(a;z x ) 



(3.15) 



W"(a;z x ) + Y,(a;z x )q 2 + ■ ■ ■ v 

and recall that a = m(z x ). 

We can easily generalise the analysis to derive ef- 
fective Hamiltonians containing any number of spin- 
and interfacial-like collective coordinates by imposing 
spin and crossing constraints at arbitrary positions. Of 
particular interest is the two-field model H[a,£] with 
cr(y) the d — 1 dimensional plane of spins at the wall 
and £(y) the collective coordinate for the surface of fixed 
magnetization m x = 0. This is a coupled Hamiltonian 
similar to the H[£\,£ 2 \ theory describing the interaction 
between order parameter fluctuations near the interface 
and wall. The model can be constructed from the 
properties of the planar constrained profile m 7r (z; a n ,£ % ) 
which satisfies the usual Euler-Lagrange equation subject 
to the constraints 



m w (0; ov, 4-) = ov 

m,r(4-;ff7r,^r) = 



(3.19) 
(3.20) 



using now standard methods. We simply quote the 
results since the details of a more involved calculation 
will be discussed later. The Hamiltonian is 

H[a,l]= J dy<jj£ u (Vcx) 2 + £ 12 V<r.W 



+ 2^22(V^) 2 + W 2 (<7j) 



(3.21) 



where explicit evaluation of the various terms in the 
reliable double parabola approximation for 4>(m) yields 



W 2 {a,t) = -r(a - a ) 2 + 2Km a (o - m a )e 



-Ki 



+k (m 2 a + (<r — m a ) 2 ) e 



2\ „-2k£ 



(3.22) 



with 



Following our earlier observations for the FJ Hamilto- 
nians let us focus on the prediction for the correlation 
exactly at the position of th e spin-plane for the model in 
question. In this case (3.14) reduces to 



^0 ' ( zX ' zX ) 



1 



W"(a;z x ) 



(3.16) 



which can be explicitly calculated using techniques very 
similar to those described in Sec. 2.2 of |l5 |. This 
procedure identically reproduces the analytic expression 
for the MF correlation function so that we identify 



G y(z x ,z x ;z x )=G (z x ,z x ) 



(3.17) 



In this way we can reconstruct the MF correlation 
function from a set of spin-like Hamiltonians {H [a; z x ]} 
using each element at the appropriate local position in 
space. Similar remarks can be shown to be valid to order 
q 2 ( and beyond if rigidity-like terms are allowed for in 
(|3.7|)) so we may identify 



r = c + n 

(To = m a + t (3.23) 
Similarly the stiffness coefficients are given by 



Eli 



dz 

o V Cov 



<y(y)Av) 



1 

2^ 



El2 



dz 

o V Cov 



m a {n£ - l)e~ Kt + 0(e 



<r(y)Ay) 



(3.24) 



(3.25) 



and 



°{y)Ay) 

-Kt 



S a/3 + 2Km a (a - m a )e Kl - 2n 2 m 2 a £e~ 2Kl 
+k [2(<j - m a f + 3m 2 a ] e" 2 ^ (3.26) 
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Note that minimization of W(a,£) recovers the MF 
surface magnetization and position of the a(3 interface. 
The two-field model H[cr,£] is an alternative candidate 
for describing coupling effects at wetting transitions and 
has its own CFRS which enable us to address problem 
(P3). We again simply quote the results since a more 
general proof will be given in the next section. 

First construct the appropriate direct correlation 
function matrix 



C(q) 



where this time 



r) 2 
°11 

fl 2 
°21 



fl 2 
°12 

f) 2 
°22 



W 2 (a,£)+q 2 -E 



(3.27) 



r) 2 — 92 ■ r) 2 — f) 2 — ° 2 ■ fl 2 — ° 2 (3 2R) 
U ll — da 2 ' U 12 — u 21 — dadl ' 22 — dp \o.^O) 

and all derivatives are evaluated at equilibrium. The 
required MF order parameter correlation functions are 
then given by 



G(0,0;q) = S u (q) 
G(0,l;q) =-m'(f)Si 2 (q) 
G(£ 1 £;q)=m'(£) 2 S 22 (c l ) 



(3.29) 



where the structure factor matrix S(q) (defined analo- 
gous to ( |2. 19 )) is the inverse of C(q) and rh{£) = 0. 

Let us finish this subsection by using the H[a,£] 
Hamiltonian to calculate the non-Lorentzian form of 
G(0,0;q) for critical wetting along route (i). For 
simplicity we shall ignore the position dependence of the 
stiffness coeffi cients although these can easily be included 
(see Sec. VB for further discussion). 

Th e eleme nts o f the matrix C can be calculated from 
( |3.22 ) and (3.23). To the appropriate order we find 



(where v 



C — K 
C+K 



C u (0) =c + k- 

G 12 (0) = 2k 2 W 
G 22 (0) = 2kVt 2 



2kv 2 t 2 



and so using (3.2G) derive 



G(0,0;q) 



r -I- k 4- -21 — 2k 



(3.30) 

(3.31) 
(3.32) 



(3.33) 



for £ 2 = 



»-f-3 • This is of the desired form consistent with 



( 1.46 ) and the exact sum-rules ( 1.47 ) quoted earlier. 

In fact the form of G(0, 0; q) at critical wetting emerges 
much more naturally from the CFRS of H[a,£] model 
compared to that of the PB model because in the 
latter we have to explicitly include the vanishing of the 
magnetization gradient at the wall (recall ( 2.20] )). In 
the H[a,£] description the singularity of m\ is implicitly 
included in the binding potential W2 (<r, £) . This feature 
works against us somewhat if we use the H[a,£] model 



to calculate G(0,0;q) at complete wetting. While 
we can derive a formally corre ct expression for this 
quantity comparison with ( 1.42| ) is not as immediate 
as with the H[£i,£ 2 ] approach which has the natural 
inclusion of explicit rh 2 terms in its CFRS. Consequently 
there is no elegant analogue of the stiffness matrix free 
energy relation in the H[a,£] formalism of complete 
wetting. Thus, while both Hamiltonians H[cr,£] and 
H[£i, £ 2 ] yield precisely equivalent expressions for the MF 
correlation functio n, wh ich formally solves the Ornstein- 
Zernike equation (1.21), the spin- and interfacial-like 
representations of order parameter fluctuations at the 
wall are more convenient for critical and complete wetting 
respectively. We shall return to this point later. 



B. Invariance of correlations within the CFRS 

In the section above we stated that the two different 
collective coordinate models H[a,£] and H[£i,£ 2 ] consti- 
tute equivalent approaches to recovering MF correlation 
functions (near the wall and af3 interface). That is, apart 
from the question of convenience for analysing particular 
transitions both choices of collective coordinate give the 
same final expression for the physical order parameter 
correlation functions G(z M , z„;q) at MF level. Here we 
give a general proof of the invariance of the G(z M , z„; q) 
with respect to different effective Hamiltonian descrip- 
tions. This section is rather technical and the reader 



may prefer to skip to Sec. Ill C where we summarise the 
main points. 

Let us denote the collective coordinates Xi(y) and 
X 2 (y) which we use to represent order parameter 
fluctuations at the wall and interface respectively. The 
proof can be trivially extended to any number of fields 
at arbitrary positions. Associated with the collective 
coordinates is a pair of surfaces described by Zi{y) 
and ^ 2 (y). These may be viewed as the locations of 
the fluctuating fields and are specified by appropriate 
functions 



Z i (y) = Z i (x 1 (y),X 2 (y)) 



(3.34) 



which will depend on the choice of collective coordinate 
system. 

The magnetization profile is doubly constrained by the 
conditions 

m(y,^(y);X 1 ,X 2 ) = M i (x 1 (y),X 2 (yfj (3.35) 

and notice that we have allowed the value of the 
magnetization at position Zi{y) to depend on one or 
indeed both of the collective coordinates. The effective 
Hamiltonian is defined in the usual way as a constrained 
fluctuation sum 



H[X U X 2 



In 



T>me 



-H LGW [m] 



C 



(3.36) 
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We will assume that it is possible to find collective 
coordinates Xi , X 2 such that the F J saddle point 
approximation is valid. Then 



H[X 1 ,X 2 }=H LGW [m s ] 



(3.37) 



wh ere ro s is the profile that minimises Hlgw subject 
to ( 3.35 ). Th us ms satisfies the usual Euler-Lagrange 
equation ( |3.4|) together with the required boundary 
conditions (|3.5| ). Following all our other treatments 
we calculate m= perturbatively in terms of the p lana r 
constrained profile which satisfies the analogue of ( |3 .35 ) 



miriZin', Xi n , X2n) — Mi (Xi^ , X 27T ) 



(3.38) 



with = ZiiXx^^Xi-n). Following Fisher, Jin and 
Parry [B4J we write 



:(y,z;X 1 ,X 2 'j = m n (z; X x (y) , X 2 (y ) 

+e(r;X 1 ,X 2 ) (3.39) 



where £(r; X\, X2) denotes the error term. By construc- 
tion £ vanishes for planar configurations 

Xi(y) = X tn Vy (3.40) 

and also at the position of the collective coordinates 

z = Zi(y) (3.41) 

for any y. 

Now let us denote the equilibrium MF values of the 
fields X\ and X 2 respectively. These may be found 
by minimisation of the Hamiltonian H[Xi, X 2 ]. Clearly 
setting Xi — Xi recovers the MF profile 



m^(z;X 1 ,X 2 ) = fri(z) 



and together with (3.35) this implies 
m{Z l )=M l {X 1 ,X 2 ) 



(3.42) 



(3.43) 



with Zi the MF position of the collective coordinate Xi. 

To elucidate the reconstruction scheme we need to find 
a relationship between the magnetization pair correlation 
functions of the LGW theory and the collective coor- 
dinate correlations of the effective Hamiltonian. Recall 
that the former are defined by 



G(ri,r 2 ) = G(z 1 ,z 2 ;y 2 - yi) 

= (m(ri)m(r 2 )) - (m(ri))(m(r 2 )) 

and satisfy the Ornstein-Zernike equation fl52[ 

8 2 H LGW [m] 



(3.44) 



rfr 



5m(ri)6m(r) 



G(r,r 2 )=<5(n-r 2 ) (3.45) 



Similarly, for the macroscopic collective coordinate 
theory 



S(yi,y2) = S(y 2 - yi) 

= (%i)% 2 ))- Wyi)>(% 2 )) 



(3.46) 



which obey |52| 



SX ll (y 1 )8X p (y) 



X\ — X\ — x 



- S pu (y,y 2 ) = 



(3.47) 



and note our implicit use of the summation convention. 

Considering the Ornstein-Zernike direct correlation 
function 



5 2 H LGW [m~\ 



(3.48) 



with an eye on (3.37), the correlation functions <S(yi, y 2 ) 
can be introduced by changing the functional derivatives 
to be with respect to Xi rather than m^. For this we 
need to use the functional chain rule 



5m s {y,z) 

'«i(y0 



dyi< 



SX 2 ( yi ) 



so that we can write 
S 2 H LGW [m s ] 



Sm s (r) SX 2 (y! 



(3.49) 



5m-E.(r)5mE,(r') 

'^i(yi) 




s 



SX 2 ( yi ) 



Sm 3 (r) JXi(yi) 
<Wi(y 2 ) S 



6ms(r) SX 2 (yi) 
, SX 2 {y 2 ) 6 



5m s (r') SX^yz) Sm s (r<) SX 2 (y 2 ) J 
xH[Xi,X 2 ] (3.50) 
Fortunately this can be written in matrix form. Defining 



x(y,r) 
<?(yi,y 2 ) 



(5m- (r) 
SX 2 (y) 
(5 ms (r) 

S 2 H[X U X 2 



8 2 H[X U X 2 ] 



SX 1 (y 1 )SX 1 (y 2 ) 8X 1 (y 1 )8X 2 (y 2 ) 

8 2 H[X U X 2 ] 8 2 H[X l ,X 2 ] 

8X 2 {y 1 )8X 1 {y 2 ) SX 2 ( yi )8X 2 (y 2 ) J 

we can rewrite ( [3.5 
5 2 H LGW [m] 



(3.51) 
(3.52) 



8m(r)Sm(r') m=ms 
J dy 1 dy 2 x T (y 1 ,r)C(y 1 ,y 2 )x(y 2 ,r') 



(3.53) 



The elements of the effective Hamiltonian direct correla- 
tion function matrix C(yi,y 2 ) are 



C M „(yi,y 2 ) 



5 2 H[X ll X 2 ] 
5X^ yi )8X„(y 2 ) 



C MV (yi-y 2 ) (3.54) 
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with derivatives are evaluated at equilibrium Xi = Xi. 

So far we have succeeded in writin g the first part 
of the Ornstein-Zernike equation ( 3.45| ) in matrix form 
involving functions of the effective theory. It is a simple 
matter to show that the magnetization correlations can 
be expressed in a similar way. Define 



w(y, r) 



Srn-, (r) 
(5m~ (r) 



(3.55) 



then 



G(r,r')= J dy 1 dy 2 w T {y 1 ,v)S{y 1 ,y 2 )w(y 2 ,v l ) (3.56) 

where S is the functional inverse of C 

J dyC w (y 1 ,y)S piy {y,y 2 ) = <5 M „S(yi - y 2 ) (3.57) 

Clearly S is the matrix of correlation functions 
Sfj,v(yi,y 2 )- To prove (3.56) notice that the Ornstein- 
Zernike equation can be written 

5(ti - r 2 ) = J drdu 1 d\i 2 dv 1 dv 2 ac T (ui,ri)C(ui, u 2 ) 

xx(u 2 , r)w T (vi, r)5(vi, v 2 )w(v 2 , r 2 ) (3.58) 



where we have used ( 3.53 ) and ( 3.56| ). Next from the 
chain rule identity 



SX^y) Sm s (r) dr _ SX^y) 



Sm s (r) SX^x) 
one observes that 



8X u (x) 



5(y-x)6^ (3.59) 



J drx(u 2 ,r)w T (vi,r) = / 



• 6Xi(u 2 ) 
5ms (r) 
<5X 2 (u 2 ) 

5ms (r) 

(5ms 00 5ms ( r ) 
Ki(vi) <5X 2 (v!) 



1 
1 



6(u 2 - vi) 



dr 

(3.60) 



Then, from (3.57) the r.h.s. of ( 3.58| ) becomes 
duia; T (ui, ri)w(ui, r 2 ) = / dui 



«i(m) sx 2 ( Ul ) 

(5m= [ri ) 5ms (ri) 



<5rra 5 (r 2 ) 
<5Xi(ui} 
i5m s (r 2 ) 
<SX 2 (ui) 

5(ri - r 2 ) 



(3.61) 



using the identity 
fdy L SXl j y) Sn " ' 



) «Xi(y) 

<5(ri - r 2 ) 



«?^ 2 (y) <5m 5 (r 2 ) 
5m H (ri) <5X 2 (y) 



(3.62) 



Equation ( |3.6l| ) is just the l.h.s. of ( p. 58; ) and 
consequently the Ornstein-Zernike relation is satisfied. 
Therefore we have shown that 



fir 



6 2 H LG w[m} 



Sm(ri)5m(r) 



G(r,r 2 ) = <J(n-r 2 ) (3.63) 



7n— ms 



If we set Xi(y) = X j the consistency condition (3.42) 
together with ( 3 . 39| ) im ply t hat the above relation 
reduces to the required (3.45). In this way the pair 
correlation function (in MF approximation) is 



G(r,r') = J dyidy 2 w T {y u r)S{y u y 2 )w{y 2l r') 



X i= Xi 

(3.64) 



To simplify (3.64) we first functionally differentiate 
( |3.39[ ) with respect to X, L 

= S(y - y')§^(z; X l{ y),X 2 (y)) + ™® 
OA^(y') dX^ oX^y') 

(3.65) 



This implies 

G(r,r')= / d Yl dy 2 



dX Xl 



■(*) 



9X 2 , 



xS'(yi,y 2 ) 
/ 58 



O 



dX ln V Z > 



<5(y'-y 2 ) 



(3.66) 



Here the r.h.s. is evaluated at the planar values X\ = Xi 
and X 2 — X 2 implying that the error term £ vanishes. 
However it is not necessary that the -M- also vanishes. 
In fact it only does so at particular points as we now 
show. 

As remarked earlier the Fisher- Jin-Parry ansatz ( 3.39| ) 
is constructed so that the error term vanishes when z = 
Zi(y) i.e. 



£(y,z i ( y y,x u x, 



o 



(3.67) 



Recalling that Zi{y) — Zi(Xi(y), X 2 (y)) functional 
differentiation with respect to ^(x) gives 



d£, 
az 



implying 

6£(y,z) 



SZjjy) 6£(y,z) 
z =Z U y) SX ^ SX »W 



(3.68) 



=2 4 (y) 



= -*( x -y)!5-(y)- &y> z ) 



--My) 



dX u 



dz 



(3.69) 



Evaluated at the MF values of the fields Xi the r.h.s. 
disappears since £ vanishes identically and hence so do 
its partial derivatives, implying 
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6£(y,z) 



z=Zi,Xi=Xi 



d£_ 

dz 



= 



(3.70) 



--Zi.Xi=Xi 



Therefore the error term in (3.66) only vanishes for 



s'e{ii,J 2 } 



(3.71) 



Multiplying out the matrices in ( |3.66 ) and choosing z and 
z' according to the above criterion one finally arrives at 



'^(^^(^(y-y') 



(3.72) 



x u x 2 



which is the central result of the generalized CFRS. It 
follows that for any two coll ective coordinates satisfying 
the local constraints ( |3.35| ) it is possible to recover 
the underlying MF correlations of the LGW theory 
for particle positions at and between the planes Zi 
(corresponding to the MF locations of the fields). 
Connection with our earlier statements can be made 
by simply taking the Fourier transform with respect 
to the transverse variable. The relation establishes 
that the physical order parameter correlation function 
is an invariant, independent of the choice of collective 
coordinates. 

Before we use the formalism to construct effective 
Hamiltonians based on collective coordinates with both 
spin- and intcrfacial-like components we illustrate how 
the general CFRS reduces to the matrix equations for 
the models considered earlier. 



and recall that the equilibrium positions £.- L = Zi 
are solutions of mf — fh(z ). Substitution into 
the invariance relation ( 3.72[ ) and taking a Fourier 
transform yields 



G(zi,zi;q) = m'(zi) 2 5 u (q) 
G(z 2 ,z 2 ;q) = m' \z 2 ) 2 5 22 (q) 
G(z 1 , z 2 ; q) = rh' (z^m' (z 2 )S 12 {q) 



(3.77) 



as quoted earlier in Sec. [I B 

(ii) The spin-like Hamiltonian 

The collective coordinates are X\ 
and correspond to positions Z\ ■■ 



0, Z 2 _ 



respectively. From the constraint equations (3. IE) 



we can calculate the ^Sf- terms which satisfy 



(0;<7 7r ,4)=0 (3.78) 



with 



dm 



-(zJCTtt,^) 



dz 
-m'(i) 



( 2* 7 ) ^7T ) 



(3.79) 



at equilibrium (£ v = £ and ov = a). The 
CFRS equations §3$} then follow di rectly from 
the general result ( 3.72 ) if we also use ( |3.13 ). 



Examples 



C. Summary 



(i) The PB model 

For this case the collective coordinates Xi(y) = 
£i(y) are interfacial-like and may be simply identi- 
fied with positions Zi{y) 



and the planar constraint is 



(3.73) 



(3.74) 



This leads to a simplification of the terms in 



( |3.72 ). Differentiation yields 



d£ j7r 



(z; £ ln , £ 27r ) =0 for z = 4 



(3.75) 



and 

3£itt 



(z;£i 7n £ 2 Tv) 



dz 



(z;£i^,£ 27r ) 



-m'(£i) for£ iT = £i (3.76) 



The general CFRS provides the necessary mathemat- 
ical framework for a complete solution to problem (P3). 
For arbitrary collective coordinates X%, X 2 with MF 
positions z\ (near the wall) and z 2 (near the a(3 interface) 
we can construct an effective Hamiltonian using a saddle 
point identification 



H[Xi,X 2 ] = min H L Gw[m] 



(3.80) 



subject to appropriate constraints. The MF order 
parameter correlati on fu nctions (satisfying the Ornstein- 
Zernike equation (1.21)) can be recovered using the 
invariant relation (noting the summation convention) 

dm din 



dX^ 



dX u , 



(3.81) 



where m. K {z]{Xi\) are the planar constrained profiles. 
The matrix S has elements 



V(q) = / dy^-y{5X^y)5X„{0)) 



(3.82) 
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corresponding to generalized structure factors. In 
practise we will only be interested in Hamiltonians 
modelling long wavelength fluctuations in which case the 
models all have the form 

H[X U X 2 ) = J dy\h:^{X 1 ,X 2 )VX^X v 

+W 2 (X U X 2 )\ (3.83) 



with generalized stiffness coefficients and binding poten- 
tial determined by m n (z; {Xi}). The fo rmul ae for these 
are given by obvious generalisations of ( 2 . 1 3| ) and ( 2.14 ) 
replacing the £j by X±. The structure factor matrix S is 
the inverse of the direct matrix 



J dy 12 e^C(y 12 ) 



with 



C^(yi 2 ) 



S 2 H[X 1 ,X 2 ] 
5X„{ yi )8X u {y 2 ) 



(3.84) 



(3.85) 



Focusing on long wavelength fluctuations the matrix 
C(q) has elements 



C M „(q) = 



d 2 W{X!,X 2 ) 
dX tl 8X v 



+ g 2 SV(A 1 ,A 2 ) (3.86) 



ignoring rigidity terms of 0(q 4 ) which can also be 
calculated if so wished. This scheme is valid for general 
collective coordinates (Xi,X 2 ) such as the PB choice 
(ii,£ 2 ) and spin-like model (a, I) discussed at length 
earlier. If our sole problem was (P3), that is to develop a 
effective Hamiltonian theory of correlations at the wall, 
the only remaining concern would be which choice of 
coordinates is most convenient for studying critical and 
complete wetting, say. In fact we have already answered 
this. At a complete wetting transition occurring for 
temperatures T ^> T\y the CFRS is best studied using 
the PB model H[£i,£ 2 ] because the surface gradient 
is large and does not exhibit any singular as h — ► — . On 
the other hand for critical wetting it is more convenient 
to use the H[o~,£\ model because fhl x contains important 
scaling behaviour and vanishes at T — TV for h = 0. It 
is therefore better to adopt a CFRS in which the fh! 1 do 
not enter explicitly but are implicitly contained in the 
binding potential. 

These remarks are strongly suggestive that we in- 
troduce a model based on a collective coordinate X\ 
describing order parameter fluctuations near the wall 
that can smoothly interpolate between spin-like (a) 
and interfacial-like {I) representations. With such a 
continuous family of possible collective coordinates we 
will be better placed to ask which choice is optimal for 
describing coupling effects beyond MF approximation. 



IV. OPTIMAL TWO-FIELD THEORY 

A. Proper collective coordinates 

First let us state that the choice of upper field X 2 is not 
in question and that fluctuations in the local position of 
the a/3 interface will be described by the interfacial-like 
collective coordinate £(y) corresponding to the surface 
of fixed magnetization m x — 0. For a given value of 
(planar) field i n let us denote the FJ profile 



(4.1) 



which will simplify our notation. Now the FJ theory does 
not account for fluctuations near the wall and it is natural 
to introduce a collective coordinate X\ which perturbs 
their profile in this area. We will denote the collective 
coordinate s and allow it to have both interfacial- and 
spin-like components (note that while the upper field 
certainly has a spin component nothing is lost by just 
including the much much larger interfacial fluctuations). 
This may be achieved by supposing that the doubly 
constrained profile is locally translated and enhanced. 
Specifically, choose an arbitrary position z\ close to the 



wall satisfying < z\ 



< 



The value of z\ will not 



appear in our final equations concerning problems (PI) 
and (P2) and in calculations we can set it equal to zero 
without loss of any information. Next we construct a 
doubly constrained (planar) profile which satisfies 



as well as 



T^7T (£-n 







(4.2) 



(4.3) 



The components £\, a are parameterised by the single 
variable s which we suitably scale to have the dimensions 
of length (and so allowing contact to be made with the 
PB theory). By construction s will be zero if both a 
and £1 vanish. A general linear parameterisation may be 
achieved using a free variable 6 such that 



ssin<5 



and 



a = Km a s cos 6 
Dividing these we obtain 



tan S — 



i\Km a 



(4.4) 



(4.5) 



(4.6) 



which establishes that the angle S controls the relative 
magnitudes of the spin- and interfacial-like components. 
In particular if we set S — the collective coordinate 
is spin-like since there is no £\ component and the 
constraint (4.2) simply enhances the local magnetization. 
For this choice the field theory reproduces the (<J,£) 
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analysis discussed earlier. On the other hand S = § 
corresponds to a purely interracial collective coordinate 
in which s denotes the position of the surface of fixed 
magnetization m FJ \z\]l) and is essentially the same as 
the PB theory. Generally the collective coordinate has 
both spin-like and intcrfacial-like components satisfying 



(4.7) 



which illustrates the metrical properties of the linear 
parameterisation (4.4, 4.5). We will refer to such 



generalized fields as proper collective coordinates. The 
geometrical meaning of the proper collective coordinate 
and coupling angle is shown in Fig. |^. 



m(z) 




FIG. 5. Detail of the planar magnetization profiles near 
the wall in scaled units. The broken curve shows m(z) 
which incorporates a local enhancement and translation of 
the FJ profile (corresponding to the solid line). The proper 
coordinate s and angle 8 are shown. 

We note that it might be possible to choose different 
metrical factors (which scale the magnetization and the 
distance) to ensure that (4.£ ) is dimensionally correct but 
the choice used here appears the most natural. 

We could now proceed to construct a family of cou- 
pled Hamiltonians {H[s,£;S}} characterised by stiffness 
coefficients T,^ u (s,£;S) and binding potential W(s,£;S). 
In fact we are also allowed to consider models in which <5 
varies with I. In this case the set of Hamiltonians whose 
properties we wish to study is denoted 



H = {H[sJ;6(e)}} 



(4.8) 



with elements distinguished by different choices of 
continuous function S(£). 



B. The optimal model 

The equivalence of each of the Hamiltonians in TL as 
models of the MF free energy, adsorption and correlation 



functions certainly does not extend to fluctuation effects 
in d = 3. As noted earlier there is no true non- 
trivial fixed point Hamiltonian for wetting at the upper 
critical dimension and it is possible that the different 
models in d = 3 predict slightly different behaviour for 
some non-universal quantities. For example H[s,£]0] or 
equivalently H [a, £] does not show any wetting parameter 
renormalization at complete wetting. To see this note 
that the argumen ts of the exponential terms in the 



binding potential (3.22) only depend on the upper field 
and substitution into the RG equations (2.27) leads to the 
same prediction for the value of lo determining 9 as the 
CW and FJ models. This contrasts with the H[s,£; 5] 
model, equivalent t o the PB model, which yields the 
different expression ( 2.32 ). In two dimensions however 
we can expect that all the models in 7i belong to the 
same universality class and flow to a fixed point under 
the action of the RG. 

These observations generate a host of questions. For 
example, why should we prefer the PB prediction for Co 
over the lo result equivalent to CW theory? It is not 
sufficient to chose the PB description of complete wetting 
simply because it appears to fit the simulation results 
better at high temperatures. More generally we must 
ask which of the models in Tt best approximates the non- 
universality of the underlying LGW model? Recall that 
our hypothesis concerning the physical origin of problems 
(PI) and (P2) is that the CW model does not account for 
the coupling of interfacial fluctuations to the relatively 
small order parameter fluctuations near the wall. With 
this our task reduces to finding the coupling angle 5 
which most accurately models small order parameter 
fluctuations in a Gaussian ap proxi mation (which is the 
basis of the RG flow equation ( 2.27 ) ) . We will argue that 
there is a unique answer to this question. 

Let us concentrate exclusively on the order parameter 
fluctuations at the wall. To this end consider a system 
with no CW-like fluctuations which may be achieved 
by fixing the location of the planar field £ n . Now ask 
what is the correlation function near the wall? Because 
the fluctuations are small a Gaussian approximation 
certainly suffices in which case we arc led to the 
prediction (using any of the models H[s,£:5}) 



G fix ( 2l ,^;q) = 



(^(z i; M.)) : 



( s A;5) + <z 2 £ii(sA;<$) 



s=S=0 



(4.9) 

which by the invariance of the CFRS is independent of 
5. We may view this as the correct MF expression for 
the correlation function in a thin-film magnet of width 
£ with fixed spins m x = at the z = £ layer (so that 
G fix (z2, Z2] q) = 0). Now focus on the value S* which 
maximises the equilibrium curvature: 



d 2 W 



d 2 W 



(s,£„;5*) for 5 ^6* (4.10) 
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with s = 0. Note that in general 5* will depend on 
i as well as T,h\,h etc. Due to the invariance of the 
G Rx (zi, z\\ q) with respect to 5 the same value must also 
maximise the stiffness matrix element: 



E u (0,^5$) <E u (0,4;O 



(4.11) 



The same argument implies that the variation of the 
order parameter with collective coordinate s is also 
largest for this value: 



ds 



for s = 
(4.12) 



Thus with this choice of coupling parameter 6 = 5*; 

(i) The Gaussian approximation is most accurate 
since the curvature is maximal. Recall that the 
fundamental assumption of the RG analysis is that 
a Gaussian approximation adequately describes the 
small fluctuations in the order parameter near the 
wall. In the effective Hamiltonian theory H[s,£;6] 
the extent of these fluctuations is controlled by 
the curvature (s , £; 6) . Ensuring that this is 
maximal avoids problems associated with the PB 
theory near T = T\y when rh'i vanishes. 

(ii) The gradient approximation (i.e. the neglect of 
rigidity-like terms) implicit in the use of a Gaussian 
Hamiltonian for the lower field is most accurate. 
Importantly by maximising En we completely 
avoid problems encountered in the PB theory when 
the lower stiffness vanishes — the value of En 
in the optimal theory is always finite. We also 
note that by maximising the stiffness we can also 
hope to have maximised the cut-off Ai for the 
lower field which will no longer show any singular 
behaviour (at least for T < Tc). Therefore we 
have succeeded in taking any physics associated 
with the decoupling of fluctuations near T = Tw 
out of the behaviour of the cut-off which had to 
be assumed in the PB theory. As we shall see 
the physics of decoupling effects associated with 
crossover from complete to critical wetting has 
a much more elegant description in the optimal 
model. 

(iii) The trace over s accesses the largest region of the 
magnetization phase space since the variation of 
m,r with s is maximal. The simple geometrical 
interpretation of the optimal angle S* can be seen 
from Fig. || where it defines the normal to the FJ 
profile (in appropriate scaled units) at the point 
z = z±. 

It is instructive to compare the optimal angle 5* 
with the worst possible choice of coupling angle which 
minimises the curvature, stiffness and variation 



In the geometrical interpretation this corresponds to 
directing s along the tangent vector to m FJ (z;£ w ) at 
z = Z\. Thus to first-order variation of s pro duces 
no change in the profile so the constraint (4.2) docs 
not generate new configurations. Moreover it is easy to 
establish that the minimum values of 7)jf (0, S) and 
En(0,^ T ;5) are both zero. Thus for this choice of 6 the 
Gaussian approximation completely breaks down. 

In this way we propose that the best Hamiltonian 
in the set TL for describing coupling between order 
parameter fluctuations near the interface and wall is 



H[s,£] = H[sJ-5*] 



(4.13) 



where 8* satisfies the novel (and equivalent) variation 
conditions ( 4.10| )— ( Ll2| ). Moreover it is straightforward 
to derive a simple expression for the optimal coupling 
angle 8* in terms of the known FJ profile. To do this 



note that the explicit formula for the variation is 
dm 



ds 



■ (zi , S, 



s=0 



= Km a cos 8 — (z\ : £^ ) sin 5 

oz 



(4.14) 

which follows from differentiation of the constraint 



equation (4.2) remembering (4.4, |4.5| ). It is trivial to 
maximise the r.h.s. and hence derive 



!)tn 



tan 8* 



dz 



-( Zi; e) 



(4.15) 



which is indeed the equation for the direction of the 
normal in Fig. |[ 



C. Derivation of the position dependence 

While the fundamental equations for the construction 
of the optimal model H[s,£] may appear to be rather 
involved and complicated, the actual derivation of the 
Hamiltonian requires only minor modifications to previ- 
ous calculations. First we find the general expression for 
the Hamiltonian H[s,£;8] with d arb itrary then specify 
the optimal coupling angle via (4.15). 

The constraint equation for the planar constrained 
profile is 



m w (s w sin 5; s x , £„) = m FJ (0; £„) + Km a s r cos S (4.16) 

where we have set z\ = for convenience. It is a simple 
matter to perform the calculations with finite z\ but our 
final results are unchanged provided we chose < kz\ < 
1. This may be viewed as the necessary mathematical 
condition which ensures that we are including the effect 
of order parameter fluctuations that are indeed near the 
wall. The planar constrained profile satisfies the Euler- 
Lagrange equation 
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d 2 m- K 
dz 2 



(z;---) = <f>'(m n (z;---)) (4.17) 



Amf J (£) 



which we need to solve. As shown by FJ [|20| the all 
important position dependence of the various quantities 
is captured by the double parabola approximation |53| ] 

£(«;..■) = /s£(m-m a ) for z < I (4.18) 



with <j)(m) remaining arbitrary for m < 0. Also note 
that position independent terms such as the interfacial 
stiffness coefficient (s) are not actually evaluated in any 
particular approximation but are left undetermined to 
be later combined into the appropriate dimensionlcss 
wetting parameter(s). The values of these are estimated 
by other means independent of the construction of the 
effective Hamiltonian. This is of course essential to the 
whole philosophy of effective Hamiltonian methods. 

The solution to (4.17) is fully specified by the following 
boundary conditions 

• For < z < s n sin S 

<9m„ 



dz 



CTO-n-i — hi 



(4.19) 



z=0 



and (4.16) 



• For sin S < z < £ v 

^7T (^7T 5 $71 : £?t ) 



(4.20) 



and (14.16| ) 
• For z > £ n 



lim / u \ 

z^oo m 7r (z;s 7r ,lT T ) = nip 



(4.21) 
(4.22) 



From the explicit solution for m 7r (z; s v , £ n ) it is 
straightforward (but tedious) to construct 



H[s,l} = y^yjiE n (Vs) 2 + £i 2 Vs.W 



where the binding potential may be decomposed 



W 2 (s,£;S) = ~r(s-s ) 2 + W(s,£;5) 



Writing 



with 



,FJ 



m a + T + Am{ W 



(4.23) 



(4.24) 



(4.25) 



2nm a Ki 2kt 2kI i 



(4.26) 



C + K C + K 

we find (working to second-order in Amf J ) 

r = n 2 (c+ n)q 2 

-k{c+k) sin<5^2(c - n)q + KTsinS^j Amf J (£) (4.27) 
«o = (nq^Amf J (£) (4.28) 
and (for h = 0) 
W(s,£) = 

2nm a (m a kscos6 + t + Amf J {£)\er K{ - l - s sin s) 
+n(m 2 a + {m aKS cosS + T )2^ e -2 K (f-,si„5) (4 2Q) 



where the useful abbreviation 



, cos S + t sin 5 



(4.30) 



has been introduced. As usual we impose W(s, £;S) = oo 
for £ < s sin 8. 

It is worthwhile pausing at this stage to emphasise the 
important features of the binding potential. First note 
that if we set s — then we automatically recover the FJ 
potential which neglects coupling effects. Secondly the 
form of the two-field potential is similar to that found 
in the PB model. In parti cular the exponential terms in 
the interaction term ( 1.29 ) only depend on the combined 
variable £ — s sin S similar to the l 2 — £\ dependence in 
H[£i,£ 2 ]- If we set 5 = J note that this variable indeed 
reduces to the difference £ 2 — £i consistent with our earlier 
statements that in this limit the proper field is interfacial- 
like. Similarly if we set S = the arguments of the 
exponential terms are independent of s just as in the 
H[a,£] model confirming that for this case the proper 
field is spin-like. 

As can be anticipated the position dependence of the 
stiffness coefficients is a little more involved. We will 
only quote the results to the appropriate order that we 
will require in our calculations. As in the PB model 
the stiffness of the lower surface is essentially position 
independent and we do not incur any problems if we 
simply approximate 



Su(s,£; S) w Sii(0,oo;<5) 
_ nq 2 



(4.31) 



Note that if we set 5 = 5 we find En as ^5- which is 
in fact identical with the PB result (which of course it 
should be). On the other hand for 6 — the proper 
field is spin-like and the stiffness does not show any 
dependence on the scali ng va riable r. It is interesting 
to note that for small r (4.31) is the exact result found 
by solving (4.17) for arbitrary potential function </>(m) 
i.e. beyond the double parabola approximation. This is 
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because in the limit of r — > the profile near the wall 
is flat and the double parabola assumption does not lead 
to any error in the calculation. 

At leading order the off-diagonal element is 



£i2(s, l\ ft) ~ K 2 m a q(£ — s sin<5)e 



— k,{£—s sin 5) 



(4.32) 



and dominates the position dependence of the stiffness 
matrix S due to the algebraic pre-factor linear in I — 
s sin 5. For 5 = 2 ( 4.32 ) reduces to the PB expressi on an d 
satisfies the stiffness matrix free energy relation (2.23). 
For general 6 the analogue of this relation is not so elegant 
due to non-universal pre-factors which detract from the 
simplicity of the PB result. 

The stiffness coefficient for the surface of fixed magne- 
tization m x = is 



£ 22 (s,£;<5) = Z a p + 2Km a (m a ns cos 5 + T)e- K{e - ssinS) 
-2m 2 a n 2 u{l - s sin 8)e~ 2R{l - s sin 5) (4.33) 

showing the FJ-like negative next-to-leading order expo- 
nential term. Indeed setting s = recovers their result 
(1.37) precisely. 

The final expression required to complete our specifi- 
cation of the optimal model is that for the temperature, 
surface field and position dependence of the optimal angle 
5*. This can be easily calculated from the known FJ 
profile and we find 



tan 5* 



T 

m a 



2c 



2ct 



(c + n)m c 



(4.34) 



which needs to be substituted into ( 4.2S )-( 4.33| ). 



V. FLUCTUATION EFFECTS IN THE OPTIMAL 
MODEL 

A. Complete wetting 

The essential observation here is that in the approach 
to a complete wetting transition (route (iii) in Fig. the 
position dependence of the coupling angle is unimportant 
and we may write 



5* = tan 1 



(5.1) 



Within the H[s,£] model it is the size of this coupling 
angle that describes the qualitative and quantitative 
nature of the coupling between the two fluctuating fields. 
First consider complete wetting transitions occurring at 
relativ ely high temperatures corresponding to r > m a . 
From (5.1) we see 



c* n 



(5.2) 



so that the proper field is interfacial-like. Thus deep 
in the complete wetting regime the optimal model is 
essentially the same as the PB model vindicating their 
earlier study. Indeed provided we keep away from the 
crossover region described below the PB model is a 
suitable starting point for the discussion of fluctuation 
effects. We can therefore anticipate on the basis of the 
optimal model that the value of the wetting parameter 
determining the critical amplitude 9 is renormalized in 
accordance with the earlier coupled theory. We shall turn 
to this shortly. 

On the other hand as we consider complete wetting 
transitions occurring closer and closer to the critical 
wetting temperature we encounter a crossover regime in 
which the size of th e co upling angle diminishes rather 
quickly. Again from (|5.l|) we see 



T 



(5.3) 



in this region (r <C m a ). Consequently near to the 
critical wetting temperature the proper field is essentially 
spin-like and the optimal model is basically the same 



as the H[a,£] Hamiltonian considered in Sec. Ill A 



Recall our remarks made earlier that in the absence 
of any interfacial or translation component (so that 
the argument of the exponential terms in the binding 
potential only depend on the £ field) there is no increment 
to the effective value of the wetting parameter Q. 
Consequently due to the rotation of the coupling angle 
as the temperature is reduced we can also anticipate that 



lim 
T— >T" 1 



(5.4) 



confirming the suggestion of PBS on the basis of the PB 
model. 

These remarks clarify the qualitative nature of the 
decoupling of fluctuations associated with the crossover 
from complete to critical wetting. The mechanism for 
this in the optimal model is far more elegant than the 
rather heuristic arguments given in Sec. [I E using the 
PB theory. Importantly we have succeeded in deriving a 
model of coupling effects which does not suffer from any 
pathological behaviour of the lower cut-off as t — > 0. To 
see this note that the stiffness coefficient of the proper 
field in the optimal model is simply 



En = £nsec 2 £ 



where 



Km: 



(5.5) 



(5.6) 



is the stiffness coefficient at the wetting temperature. 
Note that EJY is perfectly well behaved as T — » Tw i.e. 
t — > 0. Al so r ecall our remarks made earlier when we 
noted that ( |5.6| ) is the exact analytical MF result valid for 
arbitrary <f>(m) and is not limited to the double parabola 
approximation. 
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The main benefit of constructing the optimal model 
is seen by focusing on quantitative predictions for the 
scale of the wetting par ame ter renormalization in the 
crossover region (where (|5.3|) is valid). To calculate the 
critical exponents and amplitudes for complete wetting 
it is sufficient to ignore the position dependence of the 
stiffness matrix elements and write the optimal model 



H[s,£} = J dyjiEntVs) 5 



— YjO? 

2 



-W(£-ssmS*) 



(5.7) 



where W(£) is the standard binding potential for wetting. 
This can be approximated by 



W{£) =h£ + 2nm a Te 



for £ > 



(5.8) 



familiar from earlier approaches. Recall also that the 
parameters En and r are related to the transverse 
correlation length £ wa at the wall-a interface by = 
. As has been emphasised before one of the satisfying 
features of the two-field approach is that the final 
predictions for observable critical amplitudes do not 
depend on terms like En and r separately but only on 
scaled combinations which can be related to physical 
quantities for other systems (such as the Ising model). 

By writing rj = s sin S* we see that the rescaled optimal 
model H[q, £) and the PB Hamiltonian have precisely the 
same form so that we do not have to perform any new RG 
calculations. Consequently we can simply substitute the 
appro priate optimal model parameters into expression 
(2.32) derived by Boulter and Parry |l7|] . In particular 
in the crossover regime when the coupling angle is small 
we are led to the prediction 



^ + n TTlCT + 0((T/ '"" )4) 



where 



k B Tn 2 



(5.9) 



(5.10) 



is the wetting parameter for the proper field evaluated 
at the wetting transition temperature. An important 
point to note here is that the cut-off A appearing in our 
central result ( |5.9| ) is perfectly well behaved. As we shall 
show in the next section the numerical value of E^ is 
very similar to E Q /3 (at any temperature) so that we can 
reasonably anticipate that the cut-offs for both fields (s 
and £) are the same. Alternatively we may go further 
and argue that the cut-off for the proper field is the 
same as the LGW model. This may arise asymptotically 
close to the wetting temperature as the proper field is 
then a spin-like variable. In the terminology of Sec. 
Ill A the spin- like Hamiltonians {H [er; z]} can in principle 
be exact ly c alculated by evaluation of the functional 
integral (3.2) since a is always a single valued variable. 



Thus, in principle the models {H [a; z x ]} carry the same 
information as the full LGW Hamiltonian concerning 
correlation functions along the spin-plane at z — z x . 
This contrasts with interfacial variables for which the 
corresponding Hamiltonian is only appropriate for long 
wavelength fluctuations because overhangs in the field 
are ignored. Fortunately these considerations are of no 
concern away from Tq and for the most part we can safely 
assume that the cut-off for the s and £ fields is the same 
as that present in the CW and FJ models. 

Thus the optimized theory yields a precise perturbative 
expression for the surface field or temperature depen- 
dence of the renormalizcd wetting parameter determining 
the observable critical amplitude 6. It follows that the 
increment to 9 vanishes quadratically as r — * + , a 
result that cannot be derived from the PB model. The 



coefficient of the quadratic term in (5.£) is determined 
by the new wetting parameter fl the value of which we 
need to estimate. 



B. Critical wetting 

In many respects the optimal model yields results 
for critical wetting that are supportive of the CW 
and FJ theories. However the model also allows a 
more precise exp ression of the queries raised at the 
end of Sec. [I E concerning the influence of coupling 
on the behaviour of local response functions near the 
wall. First let us state that if we neglect the position 
dependence of the stiffness coefficients then a lengthy but 
straightforward application of the renormalization group 
equations described in Sec. II C reveals the same results 



for the leading order critical behaviour as the CW model. 
In particular the critical exponent vu for the correlation 
length £|| is strongly non-universal and dependent on the 
usual (un-renormalized) wetting parameter ui. We do 
not repeat this analysis here since it is not a particularly 
illuminating calculation and the result s can anticipated 
by inspection of the binding potential (4.29) with 5 = 0. 
Similarly we can anticipate that the inclusion of the 
position dependence of the stiffness coefficients will drive 
the transition weakly first-order in the same way as the 
central result of the FJ theory. Actually this calculation 
requires a slight extension of the RG theory of two- 
field models since the scheme developed by Boulter and 
Parry (within linear and non-linear treatments of the 
interaction term W {£2 — £\)) |l7,|l^l outlined earlier does 
not cater for the presence of al the position dependent 
stiffnesses. Nevertheless such a scheme has recently 
been derived within a linear RG approximation to the 
interaction term and yields a very similar flow equation 
for the effective binding potential to the FJ equations 
Ipi)] . In particular it is no surprise to learn that the 
flows of W and the AE MI , mix under RG iteration 
so that the critical wetting transition is destabilised 
due to presence of the FJ-like negative next-to-leading 
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order exponential term in the £22 stiffness coefficient. 
However as is now known the value of the tricritical 
wetting parameter to* at which the transition returns 
to being second-order is rather sensitive to the proper 
inclusion of the hard-wall contribution to the binding 
potential p2| , ^8| . Consequently a detailed study of 
fluctuation-induced first-order behaviour in the optimal 
model requires a fully non-linear RG analysis including 
the position dependence of the stiffness matrix. Such 
a scheme has not been developed yet although we 
are aware of work in progress p5[ . Consequently we 
content ourselves with the qualitative but nonc-thc- 
less important remarks that the optimal model yields 
predictions for the true asymptotic critical singularities 
that are similar to the FJ and CW models. 

However there is behaviour associated with coupling 
effects that is not included in either of these models 
independent of the predictions for true asymptotic 
criticality. This is most simply seen by calculating the 
Ginzburg criterion for the surface susceptibility x i an d 
comparing the result with that of the CW model ( 1.40 ). 
The Ginzburg criterion tests the self-consistency of the 
MF approximation and yields an estimate for the size 
of the correlation length when crossover from MF to 
non-classical behaviour occurs. Now this non-classical 
behaviour need not need be specified and may correspond 
to non-universality or indeed crossover to fluctuation- 
induced first-order effects. However we note that the 
estimates given by FJ for the size of the correlation 
length when the film jumps to infinity, characteristic of a 
first-order transition, are truly enormous and will turn 
out to be much larger than the crossover length that 
we calculate below. Consequently we feel justified in 
ignoring the position dependence of the stiffness matrix 
elements in our calculation. 

The calculation is for the susceptibility xi along the 
critical isotherm (route (ii) in Fig. |l|) and follows the 
same method as that outlined for the CW model by 
Halpin-Healey and Brezin J39|. As r is small we need 
to use the full expression for th e pos ition dependence of 
the coupling angle S* given in ( 4.34 ). The Hamiltonian 
is 



H[s,£] = |dyji^(V s ) 2 + i£ Q/3 (W) 2 



+ V(s,£) 



(5.11) 



where we have written 

V(s,l) = W 2 (s,t;5*(£)) (5.12) 

and c hang ed notation to avoid a plethora of PF's. Now 
from ( 4.34| ) we note that 



-«(/-.ain««) „ e -«(<-«f) (! + o(se- Kt )) (5.13) 



where p 2 = to 2 + r 2 . In this way (and using ( 4.26 ) and 
(4.29)) we specify the explicit position dependence of the 
optimal model binding potential by 



V(s,£) = 

-k 2 {c + n)p 2 s 2 + o(l - s -) e -^- s ?) 
2 ' p 



with 



^ + h{£- s T -)-h l 2cm \^ e 
p (c + Kjp d 



a = 2nm a T 

b = k(t 2 + um 2 a ) 



(5.14) 



(5.15) 
(5.16) 



and (as usual) v = S - T ^. Terms of order 

V ' C+K 

0(se~ 2Ke , s 2 e~ Kt ) have been ignored. 

Only a Gaussian approximation is used to include 
fluctuations and the onset of non-classical behaviour is 
heralded by the breakdo wn o f the self-consistency of this 
approach. Expanding ( |5.14 ) about its minimum value 
gives 

V {s,i)~v{s,i) + \{i-if 8 ' } 



di 2 



+ 2 ( ^ S) ^ 



with of course the MF fields given by 



dV 
~di 



dV 
ds 



= 



(5.17) 



(5.18) 



Next we write the Hamiltonian in Fourier representation 
in anticipation of the functional integrals 



H[s,£] 



d 2 q (1 



(27T) 



-(E^V + ^Mq) 



--(X aP q 2 +V U )\l(q)\ i 



(5.19) 



where we have written 



V nn with all derivatives 



evaluated at s — s and 1=1. The partition function is 

Z = J VtDs exp(-H[s,£]) (5.20) 

and is evaluated in MF and Gaussian approximation. 
The MF result is 



while in Gaussian approximation 

ft, 



(5.21) 



Z=H jvsVl exp - J A^|I(S^ 2 + ^)|5(q) 



+-(E Q/3(? 2 + V«)|%)| ! 



= N exp 



■In 



d 2 q 
W) 2 



In 



A 2 {Y l f 1 q 2 + V Si 



A 2 (E af3 q 2 + V u ) 



(5.22) 
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where AT is a suitable normalization factor which plays 
no role in the following calculation. In going between 
the last two equations we have made use of standard 
functional integrals ||(| . Therefore the effective potential 
in Gaussian theory is 



d 2 q 



2 J (2tt) 2 
+ ln(£ Q(3 g 2 + V U ) 



HY™ q 2 + Vsi 



while the MF susceptibility is defined as 
d 2 V 



Xi 



dhidh 

the analogous Gaussian result is 



(M) 



(5.23) 



(5.24) 



,eff 



1 d 2 



Xi = Xi - o 



d 2 q 

2 dhdhi J (2tt) 2 
ln(E a/3 g 2 + V u ) 



ln(EjV + V ss ) 



1 



d 2 q 



d 2 V 31 



av SB dv ES 

dh ' dhi 



dhidh 

U ~2.l (2tt) 2 I Y^q 2 + Vss ~ (S^g 2 + ^ ss ) 2 



+ 



dhidh 



dVg dVu 
dh ' dhi 



Y af3 q 2 + V u iT, aj3 q 2 + V u f 



(5.25) 



Evaluation of the partial derivatives in ( 5.25 ) is very 
algebraically intense and for our calculations we used 
the computer package MATHEMATICA. Once these are 
determined however the integrals are straightforward and 
we find 



jC-l-JlhCl+tjAVj. 



1 



^A 2 



- 1 



+ i ^ln(l + £LA 2 ) + 0(^ 4 ) 



(5.26) 



where we have used the definitions £, w 2 = ■!=■#■ an d 



Wu 



The higher order terms are tiny in the 



II " £c3- 

crossover region and can be safely ignored. Following 
Halpin-Healey and Brezin we suppose that a suitable 
definition of crossover is when 



xf 



= 1 



(5.27) 



from which follows our central result 
cjlln(l+£ 2 A 2 ) + T _^-l 
l + filn(l+eLA 2 ) 



(5.28) 



In the above expressions note that A corresponds to 
the momentum cut-off for the £ and s fields and may 



be safely identified with the cut-off for the CW and FJ 
models away from the bulk critical region. 

Importantly the optimal model expression for the 
Ginzburg criterion shows the presence of the wetting pa- 
rameter n associated with the fluctuations of the proper 
field. Note that in the corresponding CW expression 



(1.40) there is no such term (which is equivalent to setting 
ft = in ( 5.28| )) with the consequence that the true 
critical regime predicted by the present optimal coupled 
model is much smaller than previously expected. As we 
shall show in the next section the value of the second 
wetting parameter is about unity for Ising-like systems 
with the result that the crossover length £|| is ver y large. 
For this case we may further approximate ( |5.28 ) by the 
expression 



A^e^ 1 (l + A 2 CLr + 



(5.29) 



Equations (fS^), ( p.28[ ) and ( |5.29| ) are the main results 
of this paper. In the next section we turn to a discussion 
of the numerical value of the new wetting parameter 
n which will allow us to make better contact with the 
anomalous simulation data at the heart of problems (PI) 
and (P2). 



C. Value of the second wetting parameter Q 

Our final task is to specify the temperature dependence 
of the second wetting parameter Q in ord er to fully 
quantify our central predictions (5J3) and ( 5.28 ). As 
with the CW and FJ theories we should regard £1 (and 
lj) as inputs into the model which must be calculated 
by reliable means beyond MF approximation (at least in 
the bulk critical region). In particular we can identify 
two immediate concerns 



(i) Establish a link between Vi and observable thermo- 
dynamic and correlation function properties of the 
underlying microscopic model. 

(ii) Determine the critical behaviour of Q as Tyy — * Tq . 

We should emphasise here that these tasks stand 
somewhat apart from the rest of our paper and will 
probably require further work. 



1. Connection with thermodynamics 

A link between the proper field wetting parameter Vt 
and thermodynamics follows from the CFRS. To see this 
note that we can use a proper field s to describe the 
surface correlations at the wall-a interface — equivalent 
to setting £ = oo in the Hamiltonian H[s,£] (with h = 
of course). According to the CFRS 
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G^(0,0;q) 



:^(0;*,oo))' 



%^(s,oc;5) +(? 2 Eii(s,oo;<5) + 



s=0 

(5.30) 



and is valid for any value of 8. Now exactly at the wetting 
transition temperature the optimal angle 8* = (with 
I = oo) and using this we find 



G™ a (0,0;q) = 



„2yW 
1 ^11 



(5.31) 



using the fact that the surface and bulk magnetizations 
are equal at T — Tyy and h = 0. From an exact sum-rule 
identity E^j we can relate the zeroth moment Go(0, 0) to 
the wall susceptibility 



Go (0,0) 



dmi 



Xii/ksT 



(5.32) 



and also recall that and 
correlation length £ v 



moment t 
write 



determine the (second- 
Consequently we can 



k R T 



2 

wet 



(a) 
Xll 



(5.33) 



and bare in mind that quantities are defined for the wall- 
et; interface. In this way the value of the new wetting 
parameter (at a particular wetting temperature) can be 
identified as 



n = 



(a) 



2 

</.'<>: 



for T = T\ 



w 



(5.34) 



In principle all the quantities on the r.h.s. can be precisely 
determined for the full LGW model. 

The argument given above is not completely rigorous 
since it is based on reinterpretation of the CFRS 
equations which has only been established at MF level. 
Thus one may legitimately worry therefore whether 
( |5.34 ) is only appropriate away from the bulk critical 
region. In case this query is well founded we will 
content ourselves that the final identification above is 
reliable outside the bulk critical regime and try a different 
method of attack for T w < T c . 



earlier |Tl|] although the final numerical estimates are 
rather similar. 

In the bulk critical region the surface correlation l ength 
£ wa is very large and the simpler expression ( 5.2S ) 
for the Ginzburg criterion is appropriate. However in 
this result the two wetting parameters do not appear 
separately but as a ratio and it is natural to focus on this 
combination. We noted before that expression ( |5.6| ) for 
Y^Yi was exact in MF theory for arbitrary <f)(m). We can 
extend this by including an arbitrary pre-factor f2{m) 
(denoted K in pi]) as the coefficient of ^(Vm) 2 in the 
LGW Hamiltonian. With this the exact result for the 
lower stiffness of the proper field at Tw becomes 



f2(m a )m 2 a K 



(5.36) 



Now following the pioneering studies of Fisk and Widom 
[57 1 (reviewed in some detail by Rowlinson and Widom 
[58 1) we can chose and <p(m) to mimic the inclusion 
of bulk critical fluctuations in a variational (i.e. MF) 
treatment of the functional HLQw[m]. This appears, 
of course, in the saddle point identification used in 
constructing effecti ve H amiltonians. From the precise 
variational result (5.36) and the definition of u> it is 
straightforward to show that the critical ratio 



ft 

-=8K 

UJ 



w 



(5.37) 



where K is the surface tension critical amplitude defined 
in (9.82) of |||. 

The quantity K naturally arises in the Fisk- Widom 
theory of the surface tension and its value has been 
estimated using a variety of techniques. Rowlinson and 
Widom quote the predictions K = | and K w | arising 
from pure MF and the rather accurate Fisk- Widom 
theory respectively. From a fundamental point of view K 
is a nice quantity because it can be calculated using the 
celebrated e expansion in the RG analysis of the LGW 
model. Therefore directly from the known first-order 
result in e we find [Ml 



V37T 1 



e + 0(e 2 ) + 



(5.38) 



2. Universality of Q 



Here we give an argument aimed at establishing the 
universality of the parameter fl as T — > i.e. for 
wetting transitions occurring close to the bulk critical 
point. In fact the argument leads to a rather precise 
estimate for the critical value 



lim 

UC =T^T 



{n(T,-..)} 



(5.35) 



complementing the known universality of u>{Tc)- The 
reasoning we give is different to the simpler one used 



evaluated at e = 1. Thus we predict 



n n « 0.92 



(5.39) 



close to the Fisk- Widom result fir 



The error in 



this estimate is small (and in the last digit) because of 
the accuracy of u>c and K. 

The arg ument given above is different to one used 
earlier |ll| which was based on a direct analysis of (5.36) 
using Fisk- Widom theory to express f% in terms of bulk 
thermodynamic functions. This also predicts that flc 
is universal but the final expression is solely in terms 
of bulk critical amplitudes and looks rather different to 
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( |5.38 ). Interestingl y how ever the final numerical estimate 
is rather close to ( 5.39 ) although we prefer the present 
argument because it reduces Qc to previously studied 
surface quantities. 



D. Connection with simulations 



1. Critical wetting; (PI) revisited 



2. Complete wetting; (P2) revisited 



Let us denote the solutions to (1.4C) and ( 5.28) ) 
determining the respective size of the MF regime in 
CW and coupled models respectively by £§7^ and £ci- 
Then we see that the effect of coupling is to dramatically 
reduce the extent of the true asymptotic regime for the 
susceptibility \i since 



sGi 

tCW 
SGi 



:i+A 2 e 



2 ) 

•wo. I 



(5.40) 



It is easiest to assume that the exponent on the r.h.s. is 
close to the critical value ^ w | . Alternatively we may 
argue that away from Tq the exponent is determined 
by the MF value of K in which case we can estimate 
~ w | . We do not know the precise value of the surface 
correlation length ^ wa in the simulation studies but a 
sensible estimate away from the bulk critical regime is 
A£u, Q « 7r leading to 



SGi 
cw 



SGi 
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(5.41) 



depending on the choice of K. The coupled model leads 
to the same prediction for the divergence of £|| as the CW 
model and we can use the estimate of Halpin-Healey and 
Brezin igtl 



0.15 



(5.42) 



for the Ising model with Tw w 0.63Tc. Here J is the 
usual spin-spin coupling strength. In this way we predict 
crossover occurs when 



h 
J 



10" 



(5.43) 



(compared to the CW result ~ 10~ 4 ) which is much 
smaller than the values studied in the simulations. 

Thus the optimised theory of coupling effects yields 
a quantitative explanation of (PI) and allows a precise 
expression of the conjectures made earlier on the basis 
of the PB model. As we have emphasised before future 
simulation studies should focus on the measurement of 
response functions local to the a(3 interface if they hope 
to see fluctuation effects consistent with the original CW 
model predictions for critical wetting. 



As regards the BLF simulations Jl 2|Jl 3| ] the most 
important results of the optimal coupled model are the 
qualitative remarks 

(i) For T ^S> Tw the PB model is recovered so that the 
effective value of the wetting parameter uj is greater 
than uj. Hence the value of 9 > 1 + 

(ii) As T — > T w the increment to uj vanishes so that 
the extrapolated value of Q (or 9) at Tyy is the 
CW result. Thus the optimise d mo del is entire! 
supportive of the prediction ( 2.48 ) of PBS 
based on the PB model. 



It is worth re-emphasising the merits of extracting the 
value of uj(Tw) by extrapolating from the data for 9 not 
least because it yields a value for uj ss 0.8 consistent 
with long standing theoretical predictions. Importantly 
it avoids the issue of whether the wetting transition is 
weakly first-order or not. Moreover the BLF data are 
taken from susceptibility measurements near to the a(3 
interface (and not at the wall) where fluctuation effects 
are very strong. The fact that the BLF simulation results 
are in quantitative agreement with these predictions is 
(currently) the most important independent confirmation 
of the hypothesis concerning coupling effects. 

Unfortunately the BLF study is not sufficiently accu- 
rate to test the quantitative prediction (5.S) although 
the data is certainly consistent with a quadratic power 
law for uj near T w (see Fig. |). The small 6* 
result simplifies further if we consider complete wetting 
transitions occurring close to the bulk critical point. If 
T is fixed close to Tc and hi > hY ', then provided that 
A£u, a and c/k are large enough we can write 



tic 



hi - h{ 



(5.44) 



which is a universal result. We accept here that forcing 
our theory to work in the vicinity of the bulk critical 
point is a contentious issue since there are well known 
difficulties interpreting effective Hamiltonian models in 
the bulk critical region (6^]. Therefore we wish to not 
dwell too long on the nature of uj near Tc- 



3. A conjecture for the value of u>c 

Before we close the paper we comment on an intriguing 
numerical coincidence that caught our attention when 
estimating the extent of the critical region for xi on the 
basis o f the coupled theory. If we focus on the Ginzburg 
result ( 5.29 ) in the bulk critical region we are naturally 
led to the evaluation of the universal number 



(5.45) 
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Taking the Fisher- Wen estimate |36|] use ~ 0.775 we find 
that this quantity is extremely close to it 2 ! It is therefore 
tempting to conjecture that the exact value of Lie is 



1 



2hi7r - 1 
= 0.7755186 



(5.46) 



although this is of course entirely speculative. 



VI. DISCUSSION 

We begin with a summary before addressing some 
pertinent questions. 

A. Summary of method and results 

The starting point of our study is the coupling 
hypothesis 

Some of the non-universal physics of wetting 
transitions at the upper critical dimension 
(for systems with short-range forces) d = 3 
is sensitive to the coupling between order 
parameter fluctuations near the a(3 interface 
and wall. 

The rest of our analysis is a thorough investigation of the 
procedures we must follow in order to construct a two- 
field model of coupling effects at critical and complete 
wetting transitions. We first address (P3) related to 
the reconstruction of MF correlations and establish the 
following result: 

For arbitrary choices of local collective coor- 
dinates X /1 (y), the two-field model H[Xi,X 2 ] 
(defined using a FJ-like saddle point identi- 
fication) can precisely recover the order pa- 
rameter correlations of the underlying LGW 
theory (in MF approximation) using the 
invariant relation 



,9m ^(z M ;{X l })|^(z I ,;{X J })5^^(q) 



dX v , 



(6.1) 



where z\ , z 2 are the MF locations of the two 
fields Xi , X 2 and S is the structure factor 
matrix. 

To consider fluctuation effects beyond MF we must 
carefully chose the collective coordinates X M (y). Fol- 
lowing all earlier studies we adopt a standard interfacial 
variable X 2 = £ for the upper field which we define as 
the surface of fixed magnetization m x — 0. The choice 
of X\ is not so obvious and we therefore focus on a set 
of Hamiltonians 



H = {H[s,£; 5}} 



(6.2) 



characterised by different coupling angles S which de- 
scribe the relative importance of the spin-like (a) and 
interfacial- like (^i) components of a proper field s with 
magnitude satisfying 



= l\ + 



(6.3) 



At the upper critical dimension the Hamiltonians in 
the set Ti can exhibit different non-universal critical 
behaviour (for some quantities). We argue that the 
model, 



H[s,£] =H[s,£;5*} 



(6.4) 



with 6* the optimal coupling angle, best describes the 
fluctuations of the order parameter near the wall as 
determined by a novel variational principle. Analysis of 
the optimal model yields the following predictions 

(1) Deep in the complete wetting regime the proper 
field is intcrfacial-likc and the model reduces to 
the PB theory. RG analysis predicts the same 
wetting parameter renormalization effect at high 
temperatures and an elegant crossover mechanism 
associated with the rotation of S* as the tempera- 
ture is reduced to Tw- This leads to the explicit 
perturbation expansion result 



u> = UJ + Q 



(T/m a ) 2 



1 + (A& 



(6.5) 



in the vicinity of Ty\ 



(2) At critical wetting the effect of coupling is to 
dramatically reduce the size of the true asymptotic 
critical regime for the local susceptibility xi a ^ 
the wall. Specifically crossover to non-classical 
behaviour should occur when the correlation length 
is 



a + A 2 e 



2 ) 



(6.6) 



much larger than corresponding CW result. 



(3) The new physics that emerges is controlled by a 
second wetting parameter Q that has no counter- 
part in simpler one-field theories. We argue that CI 
approaches a universal value Qc ~ 0.92 in the bulk 
critical regime. 

(4) The values of the critical exponents in the true 
asymptotic critical regimes of critical and complete 
wetting are the same as those of the CW theory. 
Similarly it is highly likely that the optimal model 
exhibits fluctuation-induced first-order behaviour 
for sufficiently small u> < lo* similar to the FJ 
theory. 
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B. Remaining questions 

While one can argue that the coupling hypothesis is 
justified a posteriori it is still perhaps surprising that the 
inclusion of coupling effects leads to the modification of 
the CW model results in (6.5). To finish therefore we 
try to stand back from the mathematics and give more 
physical answers to some important questions. We begin 
with 



1. Why do coupling effects appear to matter so much? 

Firstly we should emphasise that most of the physics 
is insensitive to c ouplin g effects. As mentioned in the 
last remark of Sec. VIA even though we are at the upper 
critical dimension the values of the critical exponents for 
both critical and complete wetting are identical to those 
of the CW model. Nevertheless one can force the question 
and ask why does any of the underlying physics depend 
on the coupling to a weakly fluctuating field? To answer 
we draw on recent work aimed at studying wetting at 
rough walls |n| . 

If we fix the position of the lower field ^i(y) = 
z w{y) (say) in the PB model the two-field Hamiltonian 
H[zwA2\ describes the problem of wetting at a fixed 
corrugated non-planar wall. However even at MF level it 
is now known that the surface phase diagram is effected 
strongly by minor deviations of the wall from the plane. 
In particular analytical and numerical studies of a LGW 
model show that a (planar) critical wetting transition is 
driven first-order (and occurs at a lower temperature) 
provided the width a of the walls deviation satisfies pl[| 



a > 



C + K 



(6.7) 



which is less than a bulk correlation length. More 
generally numerical studies [B2j of CW-like models of 
wetting at rough walls in d = 2 shows that the order 
of the wetting transition is altered if 



Cs>C 



(6.8) 



where £ is the roughness exponent for the free a(3 
interface and (s is the roughness exponent for the wall 
configurations. Now let us return to the case of wetting at 
planar walls but include the coupling of order parameter 
fluctuations using a two-field model. It is natural to make 
analogy between the non-planar wall and fluctuations 
near a planar wall which increase the effective area as 
seen by the upper surface. Since the fluctuations in the 
order parameter near the wall are small they can only 
correspond to some effective roughness (s = 0. Thus 
we surmise that it is possible that fluctuation effects are 
sensitive to the coupling if 



which is only marginally satisfied in d = 3. Recall that for 
systems with long-range forces MF theory is expected to 
be valid in this dimension which suggests that coupling 
can only be of any possible importance for those with 
short-range interactions. We also remark that in d = 2 
transfer matrix studies confirm that coupling does not 
change any of the universal critical properties associated 
with wetting transitions. We do not include the details 
of the calculation here since we feel it would take us 
away from our central concerns. 

These remarks complement those made earlier JfL9[ that 
in the non- linear RG analysis of the CW model |Q(a shift 
in the origin of the binding potential (i.e. the location of 
the wall) is a marginal operation in d = 3 (and irrelevant 
for d < 3). Since such fluctuations are allowed for in 
two-field theory (by the variation of s sin 5*) it is perhaps 
not surprising that some of the non-universal physics is 
sensitive to coupling effects. 

Hopefully the above discussion helps clarify the reason 
why fluctuation effects can be sensitive to coupling at the 
upper critical dimension. If we accept their importance 
it is natural to proceed to our final question 



2. How should we visual coupled fluctuations? 

In order to address this question we first recall the 
answer to a related one — what fluctuations give rise 
to the standard binding potential W(£) in the CW 
model? In the absence of long-range wall-fluid and fluid- 
fluid forces the direct interaction between the interface 
and wall also arises due to fluctuation effects. This 
issue has been emphasised by a number of authors 
notably Gompper et. al. who point out that 

within a solid-on-solid (SOS) approximation an effective 
binding potential arises due to the presence of spike-like 
fluctuations from the interface to the wall (which can be 
seen in simulation studies). 



(3-phase 




> 



3-d 



(6.9) 



FIG. 6. Example of a spike-like fluctuation in the SOS 
model carrying information from the wall to the fluid 
interface. Such configurations are believed to give rise to the 
exponentially decaying terms in the binding potential. 

Note that these are not particularly important in d = 2 
because of the contact condition £ ~ 
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In d = 3 however where £ 3> the interface 

rarely reaches the wall and the spike-like fluctuations 
are crucial to the picture, determining (at least in part) 
the structure of W(£). Beyond the SOS approximation 
we can expect the spikes to correspond to string or 
tube-like objects of /3-like magnetization that wander 
on their way from the a(3 interface to the wall. This 
is reminiscent of the Abraham-Chayes-Chayes (ACC) 
description of bulk correlations in the d = 3 Ising model 
which is known to reproduce the expected Ornstein- 
Zernike decay of the correlation function pjj |. In this 
interpretation the leading order exponential decay of 
the binding potential can be directly related to bulk- 
like correlations (i.e. wandering SOS tubes) between 
points near the aft interface and wall J6{|. The coupled 
model describes deformations to the ACC tube arising 
from interactions with local (droplet-like) fluctuations 
which increase the effective area of the wall. This is 
consistent with the heuristic scaling theory discussed in 
Sec. [ID that coupling alters the manifestation of the 



(3-phase 



direct interaction between the fluctuating surfaces. These 
ideas are illustrated schematically in Figs. [6| ^ || and only 
constitute a minor change to the overall picture implicit 
in the CW model. 




FIG. 7. In d — 2 tube or spike-like fluctuations are no 
longer needed to carry information between the wall and the 
fluid interface as the contact condition £ ~ £j_ ensures that 
the interface itself regularly visits there. 

Moreover with this interpretation it is clear why in 
d = 2 coupling is not important because the spike-like 
fluctuations are not needed to convey information from 
the wall to the interface. Nevertheless we surmise that 
coupled models are still applicable even if they do not 
change any of the critical properties associated with wet- 
ting transitions as such theories should provide a more 
reliable description of correlation function structure. 

With this picture it is also clear why it is not necessary 
to explicitly include more fields describing fluctuations at 
other positions since such effects are implicitly included 
in the effective Hamiltonian theory when we replace the 
MF k and Y, a p by their true Ising values (see Fig. |§|). 




FIG. 8. Tube-like fluctuation from the interface to the 
wall in an Ising-like system. The tube is deformed due to 
its interaction with a local droplet-like fluctuation near the 
wall typical of fluctuations at the wall-a interface. Those 
droplet-like fluctuations away from the wall or near the 
a/3 interface only result in a renormalization of the MF 
expressions for k and E a) 3 respectively. 
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